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SYNOPSIS 


TMs thesis presents a study of nuclear magnetic 
relaxation in coupled spin systems in liquids, by steady- 
state horn onucl ear double resonance. Ihe n.m.r. spectra 
of most liquids under high resolution show a multiplet 
structure. These spectra normally exhibit little informa- 
tion on spin relaxation processes operative in the system. 
The relaxation in those systems can be studied by double 
resonance experiments, in which one of the transitions in 
the spectrimi is irradiated by a strong r.f. field while 
the remaining part of the spectrum is observed by a second 
wea}c r.f. field. Irradiation causes features in the 
spectrum which depend on relaxation. An analysis of 
these features loads to significant information about 
the mechanisms of relaxation, Doiible resonance experiments 
in this work, vj'ere performed on the coupled spin systems 
of (i) ring protons of 2, 6~dibromoaniline (AB 2 spin 
system) and (ii) of 2-bromotlii azole (AB spin system). 

These experiments were performed v/ith a view to study 
the systematicG of the method and details of the 
relaxation processes. The motivation of this study and 
a brief reviev; of earlier work is given in the introductory 
chapter' of the thesis. 



The notation and the density-matrix method of 
analysis, used in this work, are presented in Chapter III, 
The inhomogeneity of the static magnetic field is found 
to play an important role in the interpretation of the 
relaxation effects and a method of including it in the 
density-matrix calculations and of determining its 
contribution to single and double resonance transitions 
is given in Sec. IIB, The. distinctive features of 
solutions of density-matrix equations for low and high 
strengths of irradiation, which have been found to be 
useful in the analysis, are discussed in Sec. IIG, along 
with a graphical method for interpretation of the spectra 
with low strengths of irradiation. 

Various mechanisms responsible for proton 
relaxation are described in Chapter III, It is often 
convenient to use the model of an isotropic external 
random field for the description of relaxation effects. 
This model includes, the overall effect of various inter- 
actions, explicit considerations of which are rather 
involved. The conditions under which various relaxation 
mechanisms can be approximated by a random field and the 
degree of correlation between these fields at the sites 
of different spins are discussed. It is shown that 
under conditions normally used in the experiments, 
several of the mechanisms can be approximated by a random 
field. Scalar coupling interaction is also examined and 
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the conditions ■under which it can he treated only as a . 
mechanism for relaxation are discussed, 

The frequency- sweep homonuclear double resonance 
experiments were performed hy adding auxiliary apparatus, 
to the Tarian HR-100 spectrometer with appropriate modi- 
fications, The details of the experimental set-up are 
given in Chapter IV, 

In Chapter V, results of double resonance 
experiments performed on the ring proton (ABg) speotr'um 
of 2, 6-dibromoaniline (20^ solution in CCl^) are given. 

The parameters for the single resonance spectrum are 

“ ’^1 89.95 + 0.1 Hz. Double 

resonance spectra obtained for various strengths and 
frequencies of irradiation show significant relaxation 
features as well as marlced variations in the ihhomogeneity 
contribution to the various double resonance transitions. 
Two mechanism of relaxation (i) isotropic random fields 
and (ii) internal dipole-dipole interactions are considered 
for detailed analysis. The results of the analysis showed 
that the relaxation is primarily due to random fields 
ha'ving a ratio of mean square fields at ''A' (position 4) 
and^'B'* (position 5 and 5) protons as 2:1, The random • 
fields at .mvi * the two ‘b** protons are found to be 
almost completely correlated while those at proton A 
and any of the B protons practically uncorrelate.d. 



Chapter VI oorLtains the proton relaxation studies 
in the iB system, of 2~hromothiazole (neat liquid). The 
single resonance parameters are = 3.55 + 0-05 Hz and 
= 18,8 + 0,1 Hz. The scalar coupling of the 
protons with which has a fast quadrupole relaxation 
rate, gives rise to hroadenirig in the speotrum and, is an 
important source of relaxation of the protons in this 
molecule. The analysis of the double resonance spectra 
obtained for various frequencies and strengths of irra- 
diation leads to some significant infoimation on the ' 
scalar coupling parameters including the magnitudes and 
relative signs of various proton nitrogen coupling 
constants, as well as the parameters describing the other 
meohanisms of relaxation in the sample. Conventional 
relaxation time and linewidtli measurements were also 
used to facilitate the double resonance analysis. The 
relaxation is found to be primarily through scalar coupling 
and internal dipole-dipole interaction. The various 
parameters obtained by the double resonance analysis are 
for scalar coupling, 9.7 Hz, Hz with 

ss^e, sign and for dipolar 
—12 

interaction r * 7.0 x 10 sec, 

o 

The thesis has two appendixes, (i) a brief 
description and listing of a typical Fortran IV computer 
program and subroutines used in the density matrix calcula- 
tions and (ii) a reprint of an earlier work by the author on 
MR spectra of an IBX 2 spin system and solvent effects on 
proton- fluorine coupling constants. 
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I - II7gE0lO3gI0IT 

BELUATIOl-GElEEAIi 

The interaction of a nucleus of magnetic moment 
a constant magnetic field k, given by"^: 

a - -VIA , . 

causes a splitting of the (21+1) energy levels of the 
nucleus, corresponding to the eigenvalues of In a 

macroscopic samples consisting of an assembly of such 
nuclei, thermal eq,uilibriiun is reached by interactions of 
the nuclei with each other and v/ith other degrees of 
freedom in the sample, laiOTm as ^'lattice”, giving rise to 
a Bolt 2 mann distribution of population amongst these 

levels and a net bulk magnetization along the direction 
of Hq, Transitions betv/een these levels can be induced 
by a radio frequency (r.f.) field of appropriate frequency, 
applied perpendicular to The transition probabilities 

for upward and domiward transitions are equal, and since 
there is an excess of population in the lower levels, this 
results in a net absorption of energy by the spin system 
from the r.f, field. This phenomenon is known as nuclear 
magnetic resonance. The absorption of energy tends to 
equalise the population of these levels and the rate of 

j- The hamiltonians in this thesis are all expressed in 
angular frequency units. 


lA SPIl 
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aT^sorption is proportional to tlie population difference. 

On the other hand, the interaction of the spin system with 
the other degrees of freedom tends to restore the eq,uilihriimi 
population difference. She latter process is knoTm as 
the relaxation process. Relaxation is manifest in the 
process hy which a macroscopic sample dewelope^ a hulk mag- 
netization IL. along the direction of the external 

O o o 

A 

field H^k (v/here is the static susceptibility) and 

that hy which the hulk magnetization in the x-y plane 
decays to zero. 

A description of the relaxation process can he 

made hy studying the time dependence of the components of 

hulk magnetization Ll(t), the equation of motion for which 

<1^ 

in the presence of a magnetic field H is given hy Bloch ^s 
phenomenological equation [l]; 


dll. ■ /, 


/V 

D 


- S 


M -M 


o 


1 - 


(I“2) 


According to this equation, for H = H_k, if iyicbCoJXu. . M 
is in any arbitrary direction, the z-component of H appro- 
aches its equilibrium value 11^ exponentially with a time 
constant kno\m as spin-lattice (or longitudinal) 
relaxation time and the x and y components of H approach 
their equilibrium value zero exponentially v;ith a time 
constant I 2 , kno\m as spin-spin (or transverse) relaxation 
time. A determination of and Q-s functions of external 
parameters can he used to derive information on the 
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interactions responsible for relaxation. In complicated 
systems, like liquid samples giving iiroltipXo-ILllie -spectra 
TindiSr metliods like double resonance are 
used for study of relaxation processes. Ihe present work 
concerns itself with, a study of relaxation processes in 
coupled spin systems in liquid phase by steady-state double 
resonance. Before describing the double resonance techn- 
ique, the particular advantage of the method for the 
study of relaxation in coupled spin systems in liquids is 
briefly outlined in the following. 

The characteristic features and the methods of 
^ experimental study of relaxation are different for different 
situations. In solids the resonance lines are usually 
broad due to dipole-dipole interactions, leading to a 
short spin-spin relaxation time Tg. The spin-lattice 
relaxation time T^ is comparatively long. The spin-spin 
relaxation may be studied experimentally through linewidths 
and line shapes of resonances v/hile, in general, both T^ 
and ^2 can be measured by experiments in vdiich the transient 
response of the spin system to an r.f, perturbation is 
monitored (e.g. spin-echo method) to determine the relaxa- 
tion times on the basis of Eq. (l-2). 

The situation is rather different, by contrast, 
in liquids. In liquids the dipole-dipole interactions are 
averaged out by the rapid molecular motions leading to 
extremely sharp resonance lines, the widths of which are 
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primarily determined by iiie inliomogeneity in the applied 
magnetic field Furthermore ^the sharpness of the resonance 
lines in liq-uids enable the detection of small variations 
in the electronic environment (chemical shift) and interac- 
tions arising from indirect spin-spin coupling in a molecule, 
which lead to a fine structure in the high resolution n4.m.r. 
spectra. Ihe frequencies and the intensities in such 
spectra can he described hy the hamiltonian^ f or a single 
molecule^ given by [ 2 ]: 

Ho = + .1. hj 

X XvJ 

(1-3) 

^vhere v . = - y. 11/211, in X7hich y. is the gyromagnetic 

wJL mXm 

■feb 

•ratio including the chemical shift of the i- •' nucleus, 
and J- . is its spin-spin coupling constant v;ith nucleus j. 
Transient experiments like the spin-echo method can be used 
for measurement of relaxation times in liquids if the 
spectra do not exhibit multiple t structure. T/hen there 
is multiplet structure, these measurements encounter 


f Exceptions to this statement occur for the resonances 
of nuclei with spin ^ 1 , relaxing due to quadrupolar 
interactionSjor for resonances of spin nuclei coupled 
to either quadrupolar nuclei or nuclei participating 
in chemical exchange etc. 
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several difficulties. Separate and T2 measurements on 
each line in a multiplet* are not usually possible because 
the lines are closely spaced, and only an overall effect 
of a niunber of lines may be .obtained, furthermorejthe 
description of the relaxation process by single exponen- 
tials with time constants and fg* O'® Bloch eg.uations, 
is inadequate, therefore, for a study of relaxation 
processes in such cases, one needs a general theory and 
an experimental method, the results of which are dependent 
on relaxation processes, \ 7 angsness and Bloch [3] and 
Bloch [4] have given a density-matrix theoretical descri- 
ption of the phenomenon of nuclear induction. The density- 
matrix treatment is quite general and suitable for 
description of relaxation effects in a variety of systems. 
Double resonance experiments, the results of which are 
interpreted by the density -matrix theory, meet the above 
requirements for relaxation studies in coupled spin systems 
in liquids, 

IB DOUBLE RESOMCE 

The double -resonance experiments consist of 
irradiation of one of the transitions in the spectrum by 
a strong r.f, field and simultaneous observation of the 
remaining parts of the spectrum by another weak r.f. field. 
The principle of this method is based on the fact that the 
strong r.f. field competes with the processes of relaxation 
in the spin system, thereby causing a variety of effects 
viz. changes in populations and lifetimes of the energy 
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levels wMclijin turn^are reflected as clianges in intensities 
and linewidths in the double-resonance spectrum [4]. file 
CLuantitative details of these features depend on the rela- 
xation processes in the spin system. Information on the 
relaxation processes is then obtained by comparing the 
observed features with those calculated for various rela- 
xation mechanisms. 

Strong irradiation also mixes various eigenstates 
of the spin hamiltonian (Eg.. 1-3) t giving rise to new 
transitions and shifts in frequencies of single resonance 
transitions [4]. These features are implicitly included 
in a detailed analysis of double-resonance spectra even 
though they do not depend on relaxation, and can be 
completely understood by a diagonalization of the spin 
hamiltonian including the irradiating field, in a frame 
rotating at the irradjating frequency [5 (Q. 

A density -matrix description of the phenomenon 
of nuclear induction as formulated by Bloch [4] and 
Redfield [6], and given by AbragaiJi [7] and Hubbard [8], 
is capable of describing all the effects of the double 
resonance experiments, including, in particular, the 
relaxation effects. In this formalism the dynamical 
behaviour of the spin system is described by rewriting 
the general equation of motion of the density matrix, 
using second order perturbation theory for the spin-lattice 
interaction. This results in differential equations for 
the various elements of the density matrix, which include 
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the effects of relaxation through the elements of a general 
'‘relaxation matrix^"^ [6] , The various relaxation matrix 
elements conta,in, for a multilevel system, the information 
implied in Bloch’s phenomenological equation 

(1-2). The results of douhle-resonaiice experiments are 
then described through changes produced in the density 
matrix due to irradiation. This involves a consideration 
of the complete relaxation matrix. A calculation of these 
changes for various probable relaxation mechanisms is then 
compared with the observed features in a steady-state 
double -resonance experiment. This method thus avoids 
some of the theoretical ambiguities and experimental 
difficulties in and T 2 measurements, mentioned earlier. 
Itirthermore these experiments provide an example in which 
a detailed consideration of the density matrix formalism 
is required and thus form an experimental test for the 
theory. 

In the following some of the previous work on 
the study of relaxation in coupled spin systems by the 
double-resonance technique,- is briefly described in order 
to discuss the motivation and the nature of the work 
presented in this thesis, 

BaldeschwLeler [5] gave a detailed method of 
calculation of relaxation effects in double resonance 
based on the density-matrix formalism of Bloch [4] and 
Redfield [6]. This method was applied to the proton double 
resonance in the single spin system of CHCl^ and some 
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discrepancies \7ere noticed between experiment and theory. 
This method of analysis was later on re-examined by Nageswara 
Rao [9] and a similar but somewhat improved method was 
suggested. It vras shovm that the approximation originally 
proposed by Bloch [4] for systems v/ith weak relaxation, 
for the analysis of the double resonance spectra obtained 
under conditions of strong irradiation, is applicable 
even when the frequency of irradiation is very close to 
a single resonance transition. A reinvestigation of proton 
double resonance in CHGl^, on this basis, removed the 
discrepancies between the theoiy and the experiment for 
the single spin -I- case [10], 

hages\7ara Rao, Baldeschwieler and Anderson [11], 
used the above method to analyse the results of steady- 
state frequency-s\7eep double resonance experiments on a 
strongly coupled two spin system (AB) fomaed by the protons 
in 2-bromo-5-chloro thiophene, under conditions of strong 
irradiation. The validity of Bloch approximation for 
strong irradiation spectra has been verified in this 
analysis and furthermore it was possible to clearly dis- 
tinguish between different relaxation mechanisms on the 
basis of a comparison between the observed and the calcula- 
ted relative intensities in the double resonance spectra. 
The relaxation in this molecule was found to be predomi- 
nantly due to correlated external random fields. Wakefield 
and Memory [12] recently made similar studies on an AB and 
an ABX type of spin system, using Bloch approximation, 



10 


and concluded that the correlated random fields fit their 
spectra hest. Steady-state double resonance experiments 
have also been performed on an AB spin system by Matsuoka 
and Kakiuchi [13] j to obtain quantitative information on 
the contribution to relaxation for various intramolecular 
and intermolecular interactions and concluded that a single 
exponential cannot be used to describe various relaxation 
rates. 

loggle performed transient as v;ell as steady-state 
double resonance experiments on the iJB spin system of 2,3- 
dibromo thiophene [14]. The transient experiraents involved 
a periodic observation of recovery of single resonance 
intensity after the saturating r.f. field is turned off- 
He also examined the correlations betxveen external random 
fields and concluded that both internal dipole-dipole 
interaction and external random fields with partial corre- 
lation are responsible for relaxation in this molecule. 
Kuhlmann and Baldeschwieler [15] made similar studies in 
difluoroethylenes. These authors used low strengths of 
irradiation such that line splittings and shifts do not 
occur and only intensity changes in the single resonance 
spectra are observed (generalised Overhauser effect [4]) • 

The spectra v/ere analysed by Using external random field 
model for various relaxation mechanisms. It was concluded 
that spin- rotation interaction is responsible for relaxa- 
tion between states of different ''.symmetry . In this case 
information regarding the relaxation rate of fluorine 
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could "be olDtained v/iiile observing only proton resonance. 

Steady-state double resonance studies have been 
extended by Uageswara Eao and Les singer [16], to weakly 
coupled systems containing a large number of spins. It 
was shorn that if the interactions causing relaxation in 
such systems satisfy certain symmetry requirements, it is 
possible to obtain quantitative information on the relaxa- 
tion by assjiming a random field model. The double reso- 
nance experiments were performed in etliylfluoride and the 
results indicated that the relaxation of fluorine is 
approximately four times as efficient as that of protons 
in this molecule. 

Closely connected with the studies mentioned above 

are some hetronuclear double resonance experiments in the 

AX spin systems of 0 ‘^-E in chloroform [17] and that of 
12 

H -D in gaseous HD [18] . Significant relaxation infor- 
mation has been obtained by using the density -matrix 
analysis for these systems. Interesting extensions of 
the steady-state and transient double resonance experi- 
ments have been made by Gordon and Baldeschwieler [19,20] 
and Gordon [21] in which the observing and the irradiating 
r.f. fields are periodically pulsed. Double resonance 
experiments have also been performed in liquid mixtures 
( intermole cular Overhauser effect) [22]. 

These experiments clearly demonstrate that the 
double resonance experiments coupled v/ith density-matrix 
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analysis are capable of yielding worthwliile infaEfflation 
on tbe relaxation processes in different types of spin 
systems and imder a variety of experimental conditions. 

IG PEBSEIT Wmi 

Tlie work presented in this thesis, on coupled 
spin systems, is "undertaken to study the systematics of 
the double resonance method and to obtain quantitative 
information on the various mechanisms of proton relaxation. 
The scope of the present work is restricted to steady- 
state double resonance experiments. The work is based on 

the study of proton-proton double resonance spectra of 

) 

(1) 1^2 spin system formed by the ring protons of 2,6- 
dibromoaniline studied as a 20 per cent solution in OCl^ 
and (2) AB spin system formed by the protons in 2-bromo- 
thiazole s"tudied as a neat liquid. Double resonance 
spectra were obtained, in each case, by irradiating each 
one of the single resonance transitions over a wide range 
of strengths of irradiation, 

Prom the analysis of the AB 2 spin system the 
following points have emerged; it is found that an explicit 
consideration of the effect of inhomogenei"ty of the 
magnetic field is of considerable importance in the 
determination of quantitative information on relaxation 
from double resonance spectra, Ereeman and Anderson [23] 
^^have shown that the inhomogeneity contributes, in general, 
different widths to different transitions in a double 
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resonance spectrui'ii. In the present work a method has been 
suggested to include this in the density matrix analysis 
of double resonance spectra and this led to a reliable 
estimate of the relaxation parameters. This effect has 
been overlooked in earlier studies of relaxation by double 
resonance [11-14,16,18,24]. Secondly^it is shown in this 
work that it is often useful to study the double resonance 
spectra over a wide range of strengths of irradiation. 

This is due to the fact that the solutions of the density- 
matrix equations for high strengths of irradiation are 
determined by the relative values of relaxation matrix 
elements while those ,for lov; strengths of irradiation 
depend on their absolute values. This distinction is of 
considerable value in distinguishing between different 
mechanisms. 

Tvro relaxation processes (i) external isotropic 
random fields and (ii) internal dipole-dipole interactions, 
along with various probable degrees of correlation in each 
case, were considered for a detailed analysis of relaxation 
effects in the ABg spin system. The correlation constants 
for the isotropic random fields are defined to include the 
possibility of complete correlation accompanied by different 
mean square fields. The analysis led to the conclusion that 
the relaxation in the spin system is primarily due to 
external random fields with a ratio 2j1 for the mean square 
fields at the proton (position 4) to those at the " 
protons (positions 3 and 5) . The random fields at the two 
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''B*' proton sites are foimd to be highly correlated while 
those at the ^A' proton site and any of the proton sites 
are completely nncorr elated."^ 

In the AB spin system of 2-bromothi azole scalar 
coupling of protons v/ith which has a fast quadrupole 

relazation, was found to modify the process of relaxation of 
the former to a large extent. This coupling causes significant 
broadening in the proton spectrum of this molecule. The 
double resonance study yielded information on the various 
parameters describing the scalar coupling interactions viz., 
the magnitudes and relative sign of the proton-N^^ coupling 
constants. The scalar coupling interaction was also 
examined in some detail and the conditions under which it can 
be treated only as a mechanism for relaxation are discussed. 

It is found that under the conditions of the experiment the 
scalar coupling interaction contributes significantly only 
to spin-spin relaxation, so that the contact of the spin 
system with the lattice should be through other interactions. 

A detailed double resonance analysis for both high and low 
amplitudes of irradiation cases led to the conclusion that 
the internal dipole-dipole interactions along with the’ 
scalar coupling account for the spin relaxation in this 
molecule. A graphical method for interpretation of the 


The resiJlts of this study are published in a paper by 
Anil Kimiar and, B.B. Nageswara Rao, Molec, Phys., 15 , 
377 (1968). 
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spectra with^ low* amplitudes of irradiation is adopted for 
this analysis.^ 

It is often convenient to use the model of an 

isotropic random field for the analysis of relaxation 

effects in the double resonance studies. This is because 

this interaction is linear in the spin operators I„, I /and 

X y 

and the lattice variables and possesses isotropy. This 
model represents the overall affect of various intramolecular 
and intermole cular interactions explicit consideration of 
?/hich is rather involved. As part of this work, the_ condi- 
tions under which different mechanisms can be approximated 
by a random- field model are examined. It is shown that 
under conditions normally used in experiments, several of 
the mechanisms can be approximated by such a model. 

The experiments were performed a.t room temperature 
and by frequency- sweep method. The experimental arrangement 
is similar to the standarad arrangements used in literature 
[23,25] with some minor modifications arising from the 
particular apparatus used. This arrangement was set-up and 
was used along with a Vaiian HR-100 n.m.r. spectrometer. 

The plan of presentation of the thesis is as 
follows. The notation and the density-matrix method used 
in this work are given in Sec. IIA. The treatment of 
inhomogeneity and its contribution to single and double 

The results of this study form part of a paper by Anil 
Kumar, H. Rama Krishna and B*D. Nageswara Rao, Molec. Rhys.,, 
1969 (in press) . 
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resonance spectra are given in Sec, IIB. The solutions of 
the density-matrix eq^uations for different strengths of 
irradiation are discussed in Sec, IIC along with a graphical 

V I 

method for interpretation of spectra for low strengths of 
irradiation. Various relaxation mechanisms are listed in 
Sec. IIIA, and the conditions under which they can he 
approximated hy a random field are examined in Sec. IIIB. 
Chapter III also contains a discussion of the conditions 
under which scalar coupling can he treated only as a mechanism 
of relaxation. The deta-ils of the experimental arrangement 
for frequency- sweep double- resonance experiments are given 
in Chapter IV. Chapter V contains the results of double 
resonance experiments on the iB 2 spin systoDi of 2,6-dihromo- 
aniline and their analysis. louhlo resonance spectra in 
this molecule shov/ed significant relaxation features and 
marked variations in the contribution of inhomogeneity to the 
linewidths. These featiires are discussed in detail. Chapter 
VI contains the results of the double-resonance experiments 
and analysis of the relaxation effects in the AB spin system 
of 2-bromothi azole. The conclusions and results of the work 
are simmarized in the end. There are two appandixes to the 
thesis (i) a brief discussion and a listing of a typical 
computer program used in the density-matrix analysis and (ii) 
a reprint copy of an earlier work on solvent effects on 
coupling constants, by the author. 
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II - THEORY 

IIA BASIC DENSITY MATRIX EORMAIISM 
1 . Equation .of Motion 

The density-matrix method of analysis of double 
resonance spectra given recently by Baldeschwieler [5] 
and lageswara Rao [9]> based on the formalism of Bloch [4] 
and Redfield [6], (also discussed by Abragam [7])j is 
discussed briefly in this section. This presents the 
basic assiimptions involved in the, theory and the notation 
to be used in the later sections.' 

The equation of motion of the density matrix 
of the total system (including the .nuclear spin system, 
the lattice and their interaction)^ in Schrbdinger repre- 
sentation^jis given by: 


II = -i[y{(t) f ], (h-d 

where ^(t) is the spin hamiltonian which is time dependent 
because of the observing and irradiating r.f. fields, 
is the hamiltonian of the lattice and is the spin- 
lattice interaction. 


Transforming Eq. (II-l) into an interaction 
representation by defining : 

= U(t) f U“^(t) 


^'^(t) = U(t) e^^^ fl' U-l(t) 


(II-2) 
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where U(t) = U(t,0) obeys the differential equation 

(11-3) 

(II-4) 

Equation (II-4) solved in two successive approximations 
leads to: 

0 

[^-^^(t-T), J^(0)]]dT. (II-5) 

A spin density matrix c is defined hys 

o^(t) = Ir^ » (II-6) 

where 

a^(t) = U(t) c7(t) U^^Ct) (II-7) 

and Tr^^l j. represents the trace over lattice states. 
Equation (II-5) is solved with the following assimptions: 
i) The time t = 0 is so chosen that the system is in 
thermal equilibrium at that time. The lattice is assumed 
to have a large heat capacity and to remain -always in 
thermal equilibrium, ^^(0). then factorizable into spin 
and lattice parts; 


dt 


n(t) flit) 


gives; 


r <? I- 

dt 


= M (-t), y ] 
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?^{ 0 ) 


x^(o) <ji(o). 


(II-8) 


V7here ^(0) = X(0) and o^(0) = dCO) are the equilibriiun 
density -mat rices at the lattice temperature for the lattice 
and spin parts respectively. 

ii) /^£'^(t) = e^"^ Imovm as relaxation hamil- 

tonian, is a random operator uncorrelated with a(0) . Ihe 
trace over the lattice states will then give zero for the 
first term on E.Il.S. of Eq.. (II-5). 

iii) The upper limit of integ^ration in Eq. (II-5) can be 

extended to infinity if t » t , where t is the correlation 

' c c 


time for time dependence of 


f (t). While t » it 

c 


should, however, be short enough to satisfy the perturbation 
expansion of Eq. (II-5) • Eor short correlation_,a (0) in 
the second term on the E.H.S. of Eq.(lI-5) may be replaced 
by o^(t). 

iv) The time independent part of is such that 

^ 1 (high temperature approximation) . This allows 
a replacement of o (t) by (cr (t) - 0(0)) in the second 
term on the E.H.O. of Eq, (II-5), while performing the 
operation Tr^^l^ . 


Transforming back to Schrodinger representation 
and noting that in liquids_> where the lattice motion may 
by adequately described classically, the operation'of 
trace after multiplying with^(0) _, can be replaced by a 
classical ensemble average [26], denoted by ^^^one 
obtains 3 
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If = a] - 


U cr(t) - a(0)]] 


/av 


d'T , (II-9) 


where 


a(0) = cTq = i (1 




1 1 


hi 


llH'i — 
n / — 




- 9 5^0' 


(II-IO) 


h is the number of spin states and q = fe./nh!I!. 

At this point one of the following two assumption 
is usually made, depending on how short is: 

(a) ' 'Extreme narrowing This assumption requires 


m 


T 

0 0 


<< 1. If it is valid_,TJ(t-T,t) of Eq. (II-9) 
becomes unity and the spectral .densities (to be defined 
later) become independent of frequency [ref. 7j p. 280]. 

(b) Eon-viscous liquid"’’ approximation. If t is not 
short enough to satisfy (a), it is most likely that 

where the summation Ir is over all the 


k 


time-dependent parts of spin hamiltonian. (In double 

ttie^e are due to 

resonancei^ the observing and irradiating r.f. fields. 

In all the cases of present interest_^ the strengths 
of these fields are very much smaller than [ref. 7, 

p. 515]). This leads to U(t-T,t) = exp (“i 
Eq. (II-9). 

Although the exact form for ^i^('t) depends upon 
the particular interaction between the spin system and 
the lattice, it can^ in general^ be expressed as: 
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(t) = ‘(-1)'^ 

L* — j 


(II - 11 ) 


where A^^^ are spin operators and are random 

functions of lattice variables expressed in irreducible 
form. The liermiticity of requires 2 

p*^^)(t) = (-l)'l f^"''l^(t) and = (-1)'^ A^”^\ (11-12) 


This form for 


#'t) 


leads to the simplification: 




/ 7 ^* 

/ p(Q.)/'-;-'\ ■p( q.) / -i-^' 


^p(q)(t) = 6^^. (^f^^^(t) f^^^(tO 

(11-13) 

known as '‘^secular approximation'*' , which arises from the 
fact that the ''non-secular terms (q ^ q'' ) vary rapidly 
with time and give zero on averaging over the ensemble[4]. 
Using Eqs. (II-ll) - (II-I5) and the *■' non-viscous liquid'"'^ 
approximation, Eq, (II- 9 ) becomes [ 26 ]; 


f| = -i[ #0, a] 


/([ b-dT 

( 11 - 14 ) 


where 


;(t-T) = ) (-1)^ A^'I^(t) E^"^^(t-T), 


q 


(11-15) 


in which 
A^ (t) = e 




( 11 - 16 ) 
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where- V q = yH^/2tc, and in writing Eq. (11-16) the coupling 
constants and the chemical shifts have been neglected 
compared to Larmor frequencies in ''extreme- 

narrqwing'*’ limit the equation of motion of o, Eq, (11-9), 
hecomes [9]^ 


-d-f = t a-®o]]/aT 

(11-17) 

2 • honl 1 e Resonance - Rotating Coordinate Transformation 


The spin hamiltonian: 

^(•t) = 3^0+ (Ii-ia) 

consists of a time independent part defined in Eq.(l-3), 

and interactions of the spins with, a wealc observing r.f . 
field ^^j_(t) and a strong irradiating r.f. field ^ 

given by; 

^^(t) = exp(icOj^t) - 1 - exp(-iC0-,_t) (11-19) 


where 


II 


- Tt ^ ''hi ^+' 


and 


^ki = - TiV2’' 


k =1,2 


( 11 - 20 ) 


( 11 - 21 ) 


in which and Vj,.jj_ are expressed in Hz, and = H^cosZ^^ti 
-H^ singj^t 3 is the observing r.f. field and H 2 = H^ 

A A 

®os W 2 ’*' ^ ^ irradiating r.f. field. 
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TraiisforLiing to a coordinate frame rotating at the 
freq.uency (-W^k) of the irradiating field through 

the operator 

ce = exp (-liOpt Z (11-22) 

3 ^ 


and noting that: 

^2 = 3^2®“^ = I>2+ + I>2- 5 

^^(t) = exp(itO t) + D^_ exp(-ioot) ^ 

[ = 0 ! [‘^ 0 ’ = 0 5 

(11-23) 


where denotes operators in the rotating frame and 
= E(i. (11-14) gives: 


i = -i[^® + «] - . fA^(t), [ 


a - oj] 


0 

d'f 






(^ - a^) 3 


(11-24) 


(11-25) 


where 




= 27t; [2l ^ 1 ^ 2 ( 1 ) + .Z-i. Jii I(i)*l(j)] t (^2+ + D2_), 


1 <D 


(11-26) 


in which 


Ai = VQj_ + 6012/2-11;. 


(11-27) 
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! (^ - a^) is tlie relaxation term defined its matrix - 
elements in a basis in which. is diagonal if the solution 

is obiained in the rotating frame through Eq. (II-24). 

I (or - o^) will similarly be the, relaxation term defined by 
its matrix elements in a basis in which is diagonal if 

the solution is obtained in the laboratory frame through 

m 

Eq. (11-14). In general in a basis a, a' .... 


aa 


(a - 






PP" 


'aa-'pp 


r ^ 

(a 




(11-28) 


where 


R 


aa-' pp 




— X 

p'a^pa “ ^ap 4^ ‘^ya^YP 




(11-29) 


in which 


^apa-'p^ " 


•foo, 


.a 


■ OO 


<^“i 


\ 


P / /av 
( 11 - 30 ) 


dr 


In writing the spectral densities 'lQ^pcc-'p^» 
the imaginary part, which is usually small, 
[ref. 7, p. 279]t. 


in this form, 
is neglected 


Ihe hermiticity of yields [9];'^'^ 


Y Ihis term is zero for the ^^extreme-narrowing^^ assumption. 

+Y These symmetry relations are obtained for^^extreme narrowing’’^ 

tc >7 

and will hold for non-viscous liquid approximation 
if the chemical shifts and coupling constants are negle- 
cted compared to X, armor frequencies in writing the Eq.( 11-16). 
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giving: 

\a'pp^ " \'ap'p ^pp'aa^ ^ ^p'pa^a ' 


(11-31) 


( 11 - 32 ) 


and 

r \app' = 

Using Eqs. II"(ll-13)j (11-15) and (II-16) we gei for 
Eq. (II-50) under ''non-viscous liquid^ approximation: 



and for extreme -narrowing''" approximation: 


-[-en 

'^apa p 




(11-35) 


5. Solution of Equation of Motion 

Equation (11-25) is solved in two steps, first 
by determining the deviation %(i) caused by the irradiating 
r.f. field ]^f^(t) in the unperturbed spin density matrix 
o^, and then by obtaining the deviation "n(t) caused by 
the observing r.f, field on the resulting system. 

The total spin density matrix, including the effect of two 
r.f. fieldSjis then expressed as; 
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= % + XCt) + 11 (-fc) ^ (11-36) 

In writing this expression, linear superposition of the 
effects of the irradiating and the observing fields on the 
density matrix is assumed. 

In the rotating coordinate system: 

c(t) = or^ + X + ^(■b) , (11-37) 


where X independent of time, since '^2 independent 
of time in the rotating frame. 


(a) Effect of Irradiation:- Setting 


1ii€ 




uy - 0 } T](t)=0, 

in the first step of^solution, Eq, (11-25) becomes in the 
steady-state: 



Oq + ^ ] - 1 



( 11 - 38 ) 


Using Eqs, (II-IO) and (II-26) and neglecting the chemical 
shift between irrad.iated nuclei compared to their JLarmor 
frequency one gets on talcing the matrix elements of Eq. (11-38) 
in the a,a^.... basis ( such that 3 ¥;i-> !a> ): 

^ -(Suiqv^) (icvg) 

4 [4(0) -4(3)]aa< - (11-39) 

Since the density matrix is hermitian E.e(Xn'fY^) “ E-eC/C /^) 

, OCOC OC Cv 

and Im ('X„„/) = -Im and the diagonal elements of 

CX Uu \X LX 
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% are real* Separating the real and imaginary parts of 
Eq.. (11-39) gives for a / a^: 




-(STiqv^) (TCV2) [1^(3) “ (11-40) 

D ' 




■*■ ^ ^aa-'pp- ^ 


0 , 


(11-41) 


and for a = a'^; 


2 


ph,H^app' 


Ee{Xpp') 



(11-42) 


Equations (II-40) to (11-42) form a set of simultaneous 
equations in the matrix elements of ’X_^of which the E 
equations given hy Eq. (11-42) are linearly dependent. 
Using the condition Ir o = Tr = 1, gives; 



(11-43) 


and a unique solution for the complete X matrix can he 
obtained. hy replacing one of the equation in (11-42) by 
Eq. (11-43). Eurther simplification of the set of Eqs. 

( 11-40) to (11-43) is possible depending on the strength 
of irradiation. This aspect will be discussed in Sec. IIC. 
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Effect o f Observing ! Pield :- The deviation in the spin 
density matrix (a^ + X) observing r.f. field ^^(t) , 

given hy h(t), is obtained using Eq.s.. (11-25), (11-37) and 
(11-38) ast 


= -1 [^0 . ^(t)] -1 [#i(t) , + X + 


(iT(t)) 


(11-44) 


A steady-state solution for rj(t) for small values of H^(t), 
small enough not to cause any saturation effects, is obtained 
by assuming that the diagonal elements of 'n(t) are zero and 
that the off-diagonal elements have a time dependence : given by 


p / 
'aa 


(■t) 


Y+ 

aa^ 


i ^ t 
e + 




-i tA- t 


(11-45) 


Using the hermiticity for p(t), substituting (II-45) in 
Eg., (II-44), collecting terms varying as exp (ito't), 
exp (-i c^yt) and constant in time, and separating the real 
and imaginary parts give a set of linear simultaneous equa- 

I’"*-' 

tions in the matrix elements of p(t). If is not very 
close to OJ^ and for small H^(t), these equations simplify 
to give [9] : 


+ n, E^e.'pp" = 


pp' 


■[%+ > 


( 11 - 46 ) 
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p<f 


■ [DV, , Im X ] - , 

‘- 1+ ’ -'aa''^ > 


( 11 - 47 ) 


in wliicli tliempper sign is important if {u>^+ ®a’’^a''^ 
small and the lower sign if (- caJ% E^-E^z ) is small. 

Only one of the off-diagonal elements of h(t), 
corresponding to the transition being observed, is important 
at a time if there are no overlapping transitions present 
in the double resonance spectrum. One then obtains: 


Re(Y- ,)= 
^ aa ^ 


( ± , ( Re X+o^ ) ] ^^/ , Im X ] 


E^ - E^.)^ -H 


( 11 - 48 ) 


-R 




/ [ D-, + , ( Re X+a^ ) ] - ( + CO +E„-E^ , ) [ D-, , Im X ] 


aa'^aa^ ^ 1 + ' 




‘aa 


i±LO+ E„ - E„,)" + R" 


'a a' ' ' ~'aa'<xa' 


( 11 - 49 ) 


^ • Signal Inten sity 

The macroscopic expectation value of an operator 
is given by a trace of the product of the operator with 
the density matrix, and therefore the signal in a nuclear 
induction experiment is given by [5f9]: 
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S= -K 


A 

dt 


Yj_ Sr 


(ly(i) a] = - i -Jj Ir{lYl)a-I.(i)a} 


rf V" r — i lA>„t i U) t _ 

~ -1-^ -=-r. \ Y_. 2r ■{ 1^(1) e ^ cT - 1_(1) e ^ ^ f 


2l. dt 47i 1^+' 

i” 


(11-50) 


where K is a proportionality constant. The diagonal 
elements of a which are time independent contribute signals 
at U} 2 t and the off --diagonal elements at ( + 6v>''+ 

{- uy^+2 (-O 2 ) > The intensity of a transition at t^)^,in 
the absorption mode^ is given by: 


^aa' “ + 




= + K'd ^ Yi (yi))„„. sin 

i 


X [ - 


-R i , 

aa ^ aa^ 


(Re X+^q) - ( + t^+S^-E^/ ) [D^_^,Im X] 


aa 


{+ u''+ E - E -f R^ / , 

^ a a ' aa '^aa' 


] . 


(11-52) 


The second terra in the numerator of Eq., (11-52) is negligible 

if Bloch approximation holds (to be discussed in Sec IIO), 

' i ' 

and the signal / becomes a lorentzian vrith half-vridth 

i -n 

^aa-^aa^ half -height. The peak signal intensity S’” 5 is 

Uu (X 

given by: 
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“f -n 

®aa' = 


+ K( >(^^6 X'^^q) ^act^ 


aa' *- + 


0" -'aa 


•p 

aa'^aa'' 


(11-53) 


where a redefined I{;(W^,H^) includes all proportionality 
constants. The integrated intensity (integrated over 
AU) from E(i. (11-52) is obtained as; 


4“ 


aa' 




(I^(i) )q.q,/[I^( i) j 



X (linewidth)^^/ 


(11-54) 


Equation (11-54) is correct even if the second term in 
the numerator of Eq. (11-52) involving Im(X) is not 
neglected, since the corresponding integral vanishes. 

This term thus superposes an intensity distribution odd 
v/ith respect to on the lorentzian, but leaves the 

integrated intensity unchanged. 

The expression for peak-height ,Eq. (II-53) , 
contains the relaxation width in the denomi- 

nator, but in a high resolution experiment the "lineshape 
is determined by a few other factors, such as the 
inhomogeneity of the magnetic field which contributes 
significantly to linewidths. These factors should be 
considered before making a comparison of the peak-heights 
with the observed spectra. This point was not explicitly 
considered in earlier studies of relaxation effects by 
double resonance [11-14,16,18,24]. 
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IIB MAGNETIC FIELD IimOMOGENEITY 
1. Sub-ensemDle Treatme n t of Density-matrix 

If tlie magnetic field is inliomogeneous* over the 
sample volume, the spin hamiltonian is different for 
different spin-systems in the sample and the treatment of 
the previous section should then Toe modified to include 
this. Abragam [ref. 7, p. 26] suggested a method of 
handling this problem by constructing sub-ensembles, such 
that the spin-hamiltonian is the same for the systems in 
each of them, but varies from one sub-ensemble to another. 
The number of systems in each sub-ensemble should be large 
enough to treat the problem statistically, but small 
enough, at the sane time to neglect variations in the 
spin-hamiltonian within a sub-ensemble. The density- 
matrix theories given earlier [4-9] and sketched in Sec. 
IIA here, v/ill then be applicable to each sub-ensemble and 
a final result can be obtained by appropriate averaging 
over all the sub-ensembles, ' 

For each sub-ensemble one can then define a 
density-matrix the equation of motion of which is 

given by; 

> (11-55) 

where is total hamiltonian for the sub- 

^ the 

ensemble and isAsame for all the spin-systems vrithin it. 
The time variation of is then given by: 
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= U.(t) ?^-(0) uT^(t) 5 (11-56) 

where U.(t) satisfies the differential equation (see Eq, 
II-3): 


i U.(t) . (11-57) 

An important point to note here is that if the density- 
matrix y for the entire ensemble is defined by averaging 
over the density-matrices of individual sub-ensembles, 
the transformation connecting ^(■t) and 5’(o) may not , in 
general, be unitary, even though that connecting fjCt) 
and ea.ch sub-ensemble is unitary [ref. 7» p.27]. 

However^in the present problem it is not necessary to 
define an average density matrix and the entire double 
resonance calciilation may be performed for each sub- 
ensemble by solving the equation of motion of T-, given 
by Eq, (II-55) , iising the theory sketched in Sec IIA, Ihe 
relevant macroscopic quantities, like the time derivatives 
of magnetization^ can also be computed for each sub-ensemble, 
and a final result may be obtained by averaging these quan- 
tities over all the sub-ensembles. If the magnetic field 
inhomogeneity may be assumed to satisfy the requirements 
of such a treatment then the signal intensity for the 
total ensemble is obtained by integrating the intensity 
obtained for a sub-ensemble given by Eq. (11-52), after 
multiplying it by the shape function of the inhomogeneity. 


dU.(t) 

dt™"^ 
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Assuming that thq shape function (normalized) 
for magnetic field inliomogeneity h(Aii>) is a lorentzian, 
given bys 


h(ACJ») = 


lAi* 

(AfO)2 + (1/2 


3 


( 11 - 58 ) 


where l/lg = y AH^ , is the half-width at half-intensity 
due to field inhomogeneity, multiplying it with Eq.. ( 11 - 52 ) 
and integrating over A from -ooto +oo , leads to a 
folding of two lorentzians, which gives another lorentzian 
having a half-v/idth eqiial to the sum of half-v/idths of 
each component [ref. .7, p. 50]. 


In addition to the relaxation width and the 
inhomogeneity width, there is an extra linewidth arising 
from the experimental set up, due to factors like minute 
instabilities in the audio-oscillators [ 23 ]. This width^ 
though small, is an unavoidable factor, Y/hich may be 
assumed to be the same for all the transitions and may 
be added to the linewidth, in the expression for intensitj^ 
in an adhoc manner. The peak-height of a non-overlapping 
transition in a double resonance spectrum is then obtained 
as: • 



k(u>i,h3_) 


aa" 


■f 






( 11 - 59 ) 
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v/here is the total half-width of the a 

UuvX ^ 

transition given hy: 


®'aa"'aa' ^ ^^aa'^^Inhomogeneity ^^^Instabiliiy , 


(11-60) 


Ihe peal!:-height in the single resonance spectnm 
is then given hy; 







2 

aa'^ 


(11-61) 


where a, a''..., is the basis in which 
^^aa^ consists of \a‘'aa''’ ^ ^^aa' 


^0 is diagonal and 
^ Inhomogeneity/ 


( ^^)lj 2 st ability^ 


2. Inhomogeneity Contribution to Linewidths 

The linevfidths in high resolution n.m.r. spectra 
of weahly relaxing spin systems are primarily determined 
by inhomogeneity in the magnetic field, For single 
resonance spectra this leads to the situation that all 
the transitions have approximately the same observed 
linewidths inspite of any variation in relaxation widths 
arising from the nature of the particular relaxation 
mechanism for the spins in the molecule. For double 
resonance spectra, however, the contribution of the 
magnetic field inhomogeneity can be different for 
different transitions as shown by Freeman and Anderson [23]. 
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Tliis is basically due to the fact that the process of irra- 
diation of one of the nuclei in a molecule^ with a monochro- 
matic frequency, causes frequency shifts of either sign in 
the transitions of other coupled nuclei, the magnitudes 'of 
these shifts depending on the ''' frequency offset* of 
irradiation. The frequency offset depends on the local 
field at the irradiated nucleus due to the inhomogeneity 
of H^, and these local fields at both the observed and 
the irradiated nuclei are exactly correlated [27]. In an 
inhomogeneous field the '''' frequency offset has a distribu- 
tion of values giving rise to a spread in the shifts of 
the line positions due to irradiation, which may either 
cancel or augment the effect of inhomogeneity, giving 
rise to lines which are narrower or broader compared to 
single resonance JLtyi&S . 


(a) Single Resonance :- In general^the width caused 

by an inhomogeneous field of width in a transition 

a — > a'', may be written as: 




Ihhomogeneity 



St^gcc' 

i 


( 11 - 62 ) 


I 1 ' 

where is the frequency of the a —>■ transition and 

gt.j^ the larmor frequencies of various spins in a molecule, 
T'or a weakly coupled single resonance spectrum Eq. (11-62) 
gives In the case of strongly coupled single 

resonance spectrum, the larmor frequencies of other nuclei 
enter the expression for through the chemical 
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sliift (gj- - CO.) and if all the strongly coupled nuclei 
telong to the same species, Eq. (11-62) essentially gives 
for all the transitions involving nuclei i" and is 
thus the same for all of them. 

(h) Eouhle Resonance - Weakly Coupled Spin Systems ;- Eor 

the double resonance case cO- enters the expression for 

. J 

OJrvAv/ (Eq. 11-62) also through the frequency offset A. of 

CX(X J 

the irradiated nucleus and for a weakly coupled spin system 
Eq. ( 11-62) hecomes: 


( Inhomogeneity ^ AH^) i (y^ 


+ Ii 




2 % 5 Aj 


(11-63) 


Thus the contribution of inhomogeneity to the linewidths 
of different transitions is different and the fractional 
change in this contribution (with reference to single 
resonance) is given by the slope of the curve oJ ^/2 tc 

Ci/lX 

versus, A.. Any transition for which this slope becomes 

u 

(-y./y.), the inhomogeneity contribution becomes zero and 
^ J 

the only width that remains is due to relaxation and 
instability. This situation has an interesting similarity 
with spin echoes in which the fanning out of the spins 
due to iuhomogeneity of the magnetic field is completely 
recovered. This also suggests the possibility that by 
double resonance the contribution of inhomogeneity to the 
linewidths in high resolution n.m.r. can be estimated and, 
therefore_^ an estimate of relaxation widths may be obtained, A 
quantitative consideration of this effect leads to an 
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accurate estimate of tiie relaxation parameters "by doulDle 
resonance. 


E stimate b y Do uble Reso na nce ;- Eor any two transitions 

'V, . . 

a a' and p p' in the d.ouble resonance spectrum of a 
weakly coupled, spin system the difference in the linewidths, 
from Eq.s. (11-63) and (11-60), is given by; 

= (Eaa'aa' “ 



> 

. 


1 



jh + A 



vi 

^3 

2 % ' 

1 


Iyj 

B i 

7 


(11-64) 


ZIHq may be readily estimated from this equation by choosing 
two transitions a — > a and p — ^ P having widely different 
slopes provided the relaxation widths are small compared 

to inhomogeneity contributions or that the tv/o transitions 

( f 

a — ^ a and p — > p have nearly equal relaxation vddths, in 
which case. the first term in Eq. (11-64) will be small 
compared to the second, and may neglected. The most 
convenient method is to look for, in the double resonance 
spectra, transitions having slopes opposite in sign and 
magnitudes close to unity, formally in double resonance 
spectra transitions with slopes opposite in sign are present. 
As shown belowj in weakly coupled spin systems^ there are 
always pairs of transitions whose slopes are equal in 
magnitude and opposite in sign. 

Consider a weakly coupled spin system of the type 
A^X^. The diagonalized hamiltonian in the rotating frame 
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for irradiation of X spins can "be written as [28]: 


. '^f 
Uo’ 


= 271 Ja^ I^(A) + [(A^ + J i^Ca))" + r 




(11-65) 


n 


m 


; I^(i) and I^(X) = 2_,. wMcii 


where I„(A) 

^ "■ 3=1 

i and 3 a^e respectively A and X types of spins. J is 

the spin-spin coupling constant between the A and the X 

nuclei. Transitions belonging to A spins are those between 

states |m(i), m(3))> |m(i) + 1, m(3)/^ • Consider a pair 

\, . 

of A transitions between the states m(i), iii(3)/’ in(i)-l» 


(3)^ and 


m(i), -m(3)) 
/ 


i(i) -1, -m(3)^ . Their 
transition frequencies are respectively by t 


~ C(A^ + J m(i))^ + m( 3) 

+ J (m(i)-l)/^ + m(3) f- (11-66) 

and 

J** 1 

= 27t^A^-[(A^ + J m(i))^ + m(3) + 

+ [ + J((m(i)-1) ^ + v|xP m(3) Jj.. . (11-67) 

The slopes of these transitions are given by: 

1 _ (A25-+Jm(i))m(3) [A^+J(m(i)-I)lm( 3) 

3^ [(A^4-Jm(i))^+v|^]'^ [^A^+J(m(i)-l)y + ^1^]^ 

( 11 - 68 ) 



41 


and 


1 p p' _ _ ii_ ^^cca'^ 

W A,. ’~ 2% ^ A j ^ 


(11-69) 


and are ecLual in magnitude and opposite in sign. It is 
then possible to estimate the inhomogeneity contribution 
-to the linewidths from the values of these slopes and the 
observed difference in linewidths. 


lie SOlUTIOHS FOR X 

<jC 

1. High Amplitude of Irradiation; Bloch Approximation 


Bloch proposed [4], that if, in the rotating 
frame in a representation in which the stationary part of 
the spin hamiltonian is diagonal, the difference between 
any two energy levels is large compared to the corresponding 
relaxation matrix elements such that : 








( 11 - 70 ) 


for each a / a"^, then the off '-diagonal matrix elements of the 
stationary part of the spin density-matrix ( X + 
negligible in that representation. This condition (11-70) 
v;as shown to obtain even when the irradiation frequency 
coincided with a single resonance transition [9], provided 



(11-71) 


Irradiation strengths termed as high* in this study refer 
to amplitudes of irradiation which satisfy the condition 
(11-71). 



42 


Under the approximation that ( X + ^ q ) is diagonal 
in the basis (known as Bloch approximation), the 

diagonal elements of 'X are given by (from Eq,. (11-42)): 


X ^aapp ^ ^ ^aapp'' (11-72) 

P P<P^ 

and the off-diagonal elements of X equal those of (-a^) . 

Th^ solution for X- obtained from Eq, (11-72) depends only 
on the relative values of various relaxation matrix elements 
and not on their absolute values. Thus the double resonance 
spectra for which the condition (II-71) is satisfied, give 
information only about the relative values of relaxation 
matrix elements. 

2, ^^Low'^^implitudes of Irradiation 


If the amplitude of irradiation does not satisfy 
the condition (11-71) , the entire set of Eqs, (11-40) to 
( 11-42) need be solved, in general. These solutions then 
depend on the absolute values of the relaxation matrix 
elements, . Two special cases may further be considered. 

If the strength of irradiation is smaller then the linewidths 
in the single resonance spectrum, only intensity changes and 
no splittings are observed in the double resonance spectrum 
(generalised Overhauser effect) [4], The strength of 
irradiation is temed as^'lcw* for such cases and it is 
convenient to obtain a steady-state solution for X('t) 4^ 
the laboratory frame rather than in the rotating frame. 
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Secondly, if in the double resonance spectrum line splittings 
are observed and the amplitude of irradiation is not high 
enough to satisfy the inequality (11-71) [this is an 
intermediate range of strength of irradiation when YH 2 
observed linewidths] , it is possible to obtain a solution 

r- 

for f in the rotating j-ra^c ' ' j without solving the 
complete set of Sqs, (11-41) and (11-42) by making what 
may be called a partial Bloch approximation [9]. (These two 
methods of solution will be discussed briefly in the 
following. 

(a) Laboratory Goordinate Solution: - (The solution for 
X(’fc) is obtained by assimiing that the diagonal elements of 
'X are independent of time and that the off-diagonal elements 
vary as ; 

i {A-Ui -i 

X('^)ab = ^ab ® + ^ab ® 3 ^ 

(11-73) 

where a, b is the basis in which is diagonal. 

The equation of motion of the spin density matrix to be 

solved for this case is given by Eq. (11-14). If in the 

single resonance spectrum there are no overlapping transitions 

present ( line ^separations >)> linewidths) ( rjimplc lino 

approximation [rof. 7, p. 522, 15]) , it can be assumod 

that the irradiation with low amplitudes, induces 

only one off-diagonal matrix element of % betwe.tri the 

of ‘ 

pair^states affected by the irradiation. If the frequency 
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of irradiation coincides ‘With a single resonance transition 
a — ^ b, one obtains, putting T)(t) = 0 in Eql. (II-14), for 
the matrix elements of ^(t) the following equations 
[ref. 7, p. 522, 15]; ' 

(II-74A) 

(II-74B) 

(II-74G) 

and 

^ ° J ' ■ (11-75) 

a 

where I « Im (zd^)^ - ^aj Xja)> , . d = 

C^l^2+ I ~ ^ “ (b|z'^|ay = /a|z''|b^ . 

Equations (II-74) and (11-75) constitute (lJ+1) homogeneous 

linear simultaneous equations in K unknowns ( and have 

a unique solution vdiich can be represented as K.qJ - 

where 1 / are functions . of R / typ-e of elements and 

a a a u u 

are independent of the strength of the irradiating r.f. 
field. The off-diagonal matrix element of X and the strength 
of irradiation are included in I. A solution for X leads 
to^ f or any transition a ^ --$>b''] 


\ab%^ ^ ^a " Xb^ ^ - 21 , 


b / a 


A ®bta'a'(>^b-^aO =-2I , 


of 


l___ ^ccdd ^ Ac 
d 

d ^ c 
c 7^a,b 


7d) 


0 , 
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■'2 '"at 


( 11 - 76 ) 


where T^/^/ = - S^/ ,3 = “lAahah ^nd q = Vm. 

The expres'sion for the intensity of any transition 
-$► b^can be obtained in a manner similar to that outlined 
in Sec. IIA, by assuming that the diagonal elements of Ti(t) aye 
zero and that the off-diagonal elements vary as exp (+icJ2_'t)« 
The expression for integrated intensity is then obtained 
as: 




(11-77) 


from which the pealc height can be obtained by, dividing by 
the total width of the line. 

For the purpose of the analysis of double resonance 
spectra it is useful to compute the ratio [S/(S'^-S)] 
vfhere S° and S are respectively the single and double 
resonance integrated intensities of the transition a^b'^ . 
This ratio is given by: 



[(c^ 


") / / _ 
o^b b 


xTTxT"^ 


1 


Vb ^2 


1 ®ab - ^a^b'' 

4 “ " '■ ■- 


|d^t 


^a'' b^ 


( 11 - 78 ) 
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Equation (II-78) gives a straigh-t line when the ratio 
[S/(S°-S) is plotted against l/ld^j, having a intercept 

which depends only on the relative values of various 
relaxation matrix elements and a slope which depends on the 
absolute values. Kuhlmann and Baldeschwieler [15] studied 
the behaviour of this ratio for a range of values of |d^ |, 
chosen such that line splittings are not observed in the 
double resonance spectra. 

For the above plot the absolute change in intensity 
of a transition betxTOen two different spectra (single and 
double resonance) is needed. This presents experimental 
difficulties unless there is a reference line in the 
spectrum which is not perturbed by irradiation. To avoid 
this^,-a comparison between the ratios of any two transitions 
a'^ — > b'^and e''-^d^ may be obtained, both in single and double 
resonance spectra. If R represents the ratio of the double 

t ^ 

resonance intensities of any two transitions a b and ■ 
c^— d'" divided by the ratio of their intensities in single 
resonance: 


E= 




o'd 


(11-79) 


then 


■R+1 




c^ d 






[T 

a b 


c a 


( 11 - 80 ) 



47 


This equation again gives a straight line having awintereept 
^“^^ah \^h^ “ which depends only 

on the relative values of various relaxation matrix elements 
and a slope which depends on the absolute 

values. 

If there are overlapping transitions present in 
the single resonance spectrum more than one off-diagonal ' 

elements of can become simultaneously significant and 
the procedure becomes involved (see Kuhlmann and 
Baldeschwieler [15]), 

(b) Partial Bloch Approximation :- If splittings are 
obtained in the double resonance spectrum, the rotating 
coordinate analysis has to be used. Hov/ever, if the strength 
of irradiation is not large enough- to satisfy the inequality 
( 11-71) for all transitions,^^. (ll-70t) , -will still be true for 
all pairs of states other than the two connected by irra- 
diation, if the other transitions have a separation, from 
the irradiated one, large compared to Vp. In such cases 
an approximation that ( + crQ)has off-diagonal elements 

only between those two states can be made [9]» and instead 
of solving all the equations (II-40) to (II-42) , a 
solution for can be obtained by solving only (lT+2) 
equation (B equations from Eq. (11-42) and one from each 
Bqs. (11-40 and 11-41). The solution for % will, in this 
case, be dependent on the absolute values; of relaxation 
matrix elements. In the limit of no splittings this solution 
reduces to the laboratory frame solution. 
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III - EBLAXATIOI^ IvIEGHMISMS MD GOERELAIIOHS 
IIIA INTERACIIOI HAMILIOlIiOJS AMD CORRELAIIOr GOKSMIS 


The liamiltonians for various interactions wMcE can 
contribute to relaxation of spin -g- nuclei in liquid samples 
and expressions for correlation constants are given in tbis 
section. A general expression for' tbe relaxation hamiltonian 
considered as a random function of time, has already been 
given in Eq. (Il-ll) and some properties of spin operators 
A^*^^ and lattice functions in Eqs. (11-12) and (11-13). 

1. External Isotropic ■ Random Elelds 

This mechanism generally includes various inter- 
molecular and intramolecular interactions for which explicit 
calculations are not feasible and a combined effect of all 
such interactions is treated as an isotropic random field 
at the sites of the spins of interest [29]. Ihe sources 
for these fields are external to the spin system of interest 
but not necessarily external to the molecule. The interaction 
of these fields with a spin 'i'" can be expressed as [ref. 7, 
p. 508]; 





Yi £(i) * 


9 


(III-l) 


such that the spin operators and lattice functions as defined 
in Eq. (II-ll) can be expressed as; 


(+ 1 ) _ 1 + 

A. = + I (i) 

^ ^2 



(III-2) 



These random fields are considered isotropic, which is 
introduced through the relation [ref. 7, p. 508, 9]; 


/ 

\ 


,(1) 




/av 




Y- 


Hf(t)j 


ST(t) 


'av 


Zl/I 

2 \ 


Hpt) 


\ 






where f(i) is the mean sauare field at the site of spin '* i'^'^. 

The degree of correlation between random fields 
at the sites of different spins in a molecufLe' may be 
defined as: 


hi' = 




av 


[<|pUJ(t)(^ >av >av3 


(III-5) 


where 0^^^/ is the correlation constant between nuclei i and 
i*^ and varies betv;een 0 and 1 for the limiting cases of no 
correlation and complete correlation between the fields 
respectively. This definition of is chosen here such 

that the mean square values of the random fields 
^|p(Q.)(t) 2 and ( t) ^ ^g_y need not be equal for 

partially or completely correlated cases as has been 
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considered ly earlier workers [ref. 7, p. 508, 11,14-16]. 

Extreme-narrov/ing asstrniption is generally used for 
calculations witli this nechanism, but depends upon the 
specific interaction being considered as a random field. 


2 . 


Internal Dipole- hip ole 


Interaction 


Ihe interaction hamiltonian for direct dipole-dipole 
coupling between different spins can be written as [ref. 7, 
p. 289, 15, 30] 




" X } (-1)^ (III-6) 


q=-2 


such that 


A. 


(± 2 ) 


•I 


+ + 
= 13 








4 °'‘ = t s„ + i; s+) ] 


"E ‘^E 4' E ^E ^ -^E “E-' 


(III-7) 


and 




(III~8) 


where the subscript E refers to a pair of spins 1^^ and 

rj;j- is the distance between them, and y| are the spherical 
harmonics of rank 2, with and specifying the orienta- 
tion of r^j- relative to the laboratory- fixed coordinate 
system. Transforming E'^^s to a molecule fixed coordinate 
system, the' correlation constants between any two pairs E 
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and can iDe olstained by averaging over all the molecular 
orientations, as [15,30]; 

1 




Js Yi .Yg (3 oos^'f^^f-1) 


■h ’"N "-r 


"If 


(III-9) 


where is the angle between r^j and r^j^ . Equation 


/ N 

(III-9) is obtained by using^addition theorem of spherical 
harmonics. For h = IT"^, = 0 and g^^^j are- Imovm as auto- 
correlation constants and for E ^ S^^r/ are called cross- 

correlation constants v/hich can be calculated from the 
geometry of the molecule. 


3. Scalar Goiipling 


Ihe scalar coLipling between a spin 1 = ^ nucleus 
and a spin S ^ 1 nucleus acts as a relaxation mechanism for 
spin I, if the spin states of are having rapid time 
dependence due to e.g. its fast quadrupole relaxation rate. 

m 

The relaxation hamiltonian for scalar coupling can be 
expressed as [ref. 7, p. 310]; 



(III-IO) 


1 

In this the spin ^ is treated as part of the lattice. The 
sma over ^ i is over all the spins I. The spin operators 
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are expressed as in Ilq. (III-2) and tiie lattice operators 
can be expressed as (see Eq, II-ll) : 

™ 1 i ( T 

E (t) = + S (t) j E^°'’(t) = S^(t) . (III-ll) 

lliere is a complete correlation between the 
interactions of various spins I(i) with the same spin S(t) 
in a molecule and it may be assumed that there is no 
correlation of these interactions with that to another 
spin S'^(t) in the molecule. 

Spin-Rotation Interaction 

The interaction of the nuclear magnetic moment 
with the magnetic field produced at its site by the rota- 
tion of the molecule containing the nucleus, is loiown as 
spin-rotation interaction and can cause nuclear magnetic 
relaxation. The hamiltonian for the spin rotation interac- 
tion is given by [31] > 

^’(t) = I • C(t) • J , (III-12) 

where J is the angular momentum of the molecule containing 
the spin I and C is the spin-rotation tensor. Assuming a 
spherical model of the molecule containing nucleus I, 
one obtains in the molecule fixed frame C = Oil and 

z 5. 

C ^ = C , = Q. and the interaction hamiltoniani" in the 
X y ^ . 

■j- A general expression for asymmetric rotor is given by 
Huntress [32] . 
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molecule fixed frame as [31]; 


imt)] 


(i) 


1 


+1 

'-■1 






(III-13) 


wiiere file supGrscrxp'i; ( i) indicafes molecule fixed frame 
values, and 


= I" ; =^-1.1- , 

^2 

j(o) ,jZ . jtiD . jyj 

^2 

and 

= % I C., =-0, . 


+1 “ 


(III-14) 


(III-15) 


(III-16) 


Transforming A and ^ from molecule fixed frame fo the 
laboratory coordinate by using the relation: 




(i) 


V ..(l) u(l) 


Z. » (III-17) 

q" 

and a similar relation for_J^ "tlie interaction liamiltonian 
in the laboratory frame is given by [31]: 

^ft) (-1)1 I-f-l) , (III-18) 


where 

p(-q) 


= 


q" , q^ 


(III-19) 
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The spin-rotation interaction is completely corre- 
lated for various spins in the same molecule and completely 
uncorrelated for spins in different molecules [31]. 

The spin-rota.tion constant 0 is usually smaller 
for protons compared to nuclei like and is usually not 
an effecient relaxation mechanism for protons ,[33j-34]. 

5. Anisotropic Chemical Shift 

Under the influence of an external magnetic field 
IJ, the induced precession of the electrons causes a inter- 
action of the type -Y I • A * H , where 4. is the chemical 
shielding tensor of rank 2 [ref. 7, p. 515]. The trace 
of this tensor gives rise to the chemical shift and the 
anisotropic part A.^ of this tensor is time dependent in 
liquids Toecause of molecular tumhling and causes relaxation 
in nuclear magnetic resonance. The interaction hamiltonian 
is given in ref. 7, p. 315. 

ho examples are laiovm in which this interaction 
makes significant contribution to proton relaxation, 
presumably because of the small chemical shifts of protons. 

IIIB CONDITIONS UITDEli WHICH VARIOUS INTERACTIONS CAN BE 
APPROXIMATED BY A RAITOOM EIELD 

The random field mechanism is often convenient to 

the 

use, forAdescription of relaxation effects^ in cases like 
the analysis of double resonance spectra. This mechanism 
is linear in its spin variables and possesses isotropy, 
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resulting in a considerable simplification in tlae calculation 
^^ 3 . large number of relaxation iaatriX' elements involved in 
tile analyses. This meciianism; serves as a model for 
describing the effects of a complex set of intermole cular and 
int^ramolecular interactions often present in a complicated 
spin system. A knowledge of the values of the random fields 
at the sites of various spins in a molecule and the correla- 
tions between these fields could then be used to indicate the 
various interactions responsible for these fields, It^ 
therefore, becomes important to examine the conditions under 
which various interactions can be approximated by a model of 
external isotropic random fields. 


General Properties of Random Pields 


The hamiltonian for external random fields and "the 
properties of their spin and lattice operators, including 
isotropy (Eq. III-4) and the secular approximation (Eq. II-I3) 
have been given earlier. The processes responsible for time 
dependence in the random field hamiltonian will be assumed to 
be Markoff ian, leading to a correlation function G(t) as 
[ref. 7, p. 297]: 




(III-20) 


where Vc is the correlation time and is the characteristic 
time for molecular reorientation, which is responsible for 
time dependence of 
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Substituting tlie form of "''^( t) from Eq., (Il-ll) , 
for a single spin, into the expression for spectral density 
of the “non-viscous liquid^^ approximation Eg, (11-34)? and 
using Eqs. (III-2) - (III-4) and the relations; 


e ° e ° e ^ 


(III-21) 


* i 

where , I = I and 1° = gives: 


+1 


^nv 


q=-l 


■ 2t / 

2 / 

[l+CcDj 


a 




(-1)'^ p(l)(t) 


av 


= f 


+1 

I 


2t 


<iT^)q TliHTT X'' 




a 


(III-22) 


In the extreme narrowing case when ^*>j% 1, (see 

X c 

Eg. 11-35): 


-j-en 


+1 

= i: 

g=-l 




2fT^[|(l )^p (I (I )^p 


+ 2(1 )q.p (I )p'^cc''^ 


(III-23) 
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Equations (III-22) and (III-23) for the spectral densities 
contain the essential features of this mechanism. The 
expressions for T^ and T2 can simply he obtained from these 
equations hy calculating the valties for a single spin ^ case 

and noting that 1/T^ = 2 Jq212 IA 2 ~ '^1122 '^1212^ 

where I^jl^ = ~ “"S'! 2^ • 

2* Various Interacti ons 

The procedure for obtaining the conditions under 
which various interactions can be a,pproximated by an iso- 
tropic random field will be to obtain the spectral densities 
under '"non- viscous liquid^ and 'extreme-narrowing ^ approxi- 
mations and compare them with Sqs. (III-22.and 23). All the 
interactions mentioned in Sec. IIIA will be considered except 
anisotropic chemical shift, which obviously cannot be approxi- 
mated by an isotropic field. 

(s.) D ipole-Eipole Int er acti on;- This is a direct magnetic 
interaction between two spins and depends both on the distance 
between the spins and the orientation of intemuclear vector , 
with respect to a space, fixed frame. The tv;o spins may either 
belong to the same molecule((i) intramolecular dipole- 
dipole interaction) or different molecule s ( (ii) inter- 
molecular dipole-dipole interaction). The time dependence 
in these dipolar interactions a,re due to molecular rotation 
in (i) and due to both translation and rotation in (ii). We 
will consider both these interactions separately. 



(i) Intramolecular dipole-dipole interaction ;- If both the 
spins of interest are spin ^ nuclei then this interaction 
is generally treated explicitly in relaxation calculations 
through what are loiown as internal dipole-dipole interactions 
the experessions for which have been given in Sec. IIIA, and 
these interactions cannot be approximated as random fields. 
However, if one of the spins has rapid time dependence in its 
spins states arising from a fast quadrupole relaxation rate 
(spin ^ 1) then this interaction, as shown belo?/,can be 
approximated as an external random field. 


The dipole-dipole interaction hamiltonian Eq.(lII- 
6-8) can be re?;ritten for a single pair as, 


/ 


f (t) 


I+[3+ Tii(~2) gS p(-l) _ ijs g- p(o)^ 


+J~ - S" + l"[-i^ ] 

+1 +1 

= ^ ^ qI gP Tp-(p+g) ^ (III-24) 

q=-l p=-l 


where. 
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-1 

-1 


= 0 °-, 
o +1 


,-l 


-1 ~ "^+1 
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0“^ = -C° = 1 

O -1 


c° = /I 

o \J 3 


(III-24A) 


Substituting this hamiltonian in the expression 
for spectral densities (Eq_. 11-34) for *^non-viscous liq,uid^ 
approximation and making use of the relations (III-21) and 
similar relations for S, one gets; 



60 


^ n n / 


^ j a^'y 


(1 a^ 

— C£:> ^ ^ / 

VfV 


KtJj q^+Wg P")V,-(pH.<i) 


4-Cp+a)(^-) p-(p'+4')(i-.^)\ . 

\ /^ 


(III-25) 


Consider now tiie case when one of the spins, say 
spin has a rapid time dependence arising from interactions 
which do not involve I* liiis can arise from rapid quadrupolar 
relaxation of spin S, if S ^ 1, The motion of spin S may 
then he assumed to be unaffected by its interaction with spin 


I and 


the rapid relaxation of S may be assumed to be 


described by Bloch equeabions. This condition is usually 
obtained for S = 1 and for S>1 if "'extreme narrowing*' is 
valid for molecular motions responsible for reltaxation of S 
[ 5]- This latter condition also leads to = t.,, the 
characteristic times in Bloch equations. In addition if 
and are very much different from t , the correlation time 
for molecular rotation, it is permissible to separately 
average the B^aiid S_£iii Eq. (I II -25) giving for St7]'. 




/s^(o) s_C) ) 

\ 4. • / 


2s(s+i: 


HtI/t- 


-|v1/to +i6L>„T 

^ G ^ f 

® * (III-26) 


Equation (III-25) v>^ill now have a correlation constant gi 


iven 
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+ 


1 


( fPl +1) 


(111-27) 


In the limit <5C j p j ^ 2_) » for most of the cases, 


-V' 


Making use of this and of Bqs. (111-20,26,27) 


and ( 11-13) , and substituting the value of from Eq, (III- 

P 

24 a) in Eq. (III-25) one obtains for spectral densities: 
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2S(S+1) 
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(III-28) 


which is identical to the expression (III-22) with 


. _ 2S(S+1) 
3 


• 


(III-29) 


Thus. the conditionsfor intramolecular dipole-dipole 
interaction to be approximated as random field are that,, the 
time dependences of spin S and that of E are widely 

-V/ 

different and that Bloch equations can be assumed for motion 
of spin S. 


Correlations : - If there are two or more spins' I (i) in the 
molecule having dipolar interaction with a spin S, then the 
correlation between the random fields produced by the 

S spin at the different spins I(i) are given by, using 
Eq. (III-9) and Eq. (III-5); 



S2 


Cii/ = i (3 - 1) (III-30) 

where '\|/'jj--/ is iiie aaigle between the vectors joining i and ±' 
spins with the spin S, In such a case it is then possible 
that the random fields nay be completely correlated (if the 
three spins lie in a straight line) having different mean 
square fields. The correlation with more than one spin 
having fast time dependence will a.lso be given by Eq. (Ill- 
30), since 

(ii) Inter-molecular dipole-dipole interaction ;- In this 
case both rotational and translational motion of the molecule 
are responsible for changes in orientation and magnitude of 
internuclear ascis. However, assuming in the following that 
the molecules are spherical end that the spins of interest 
are at centers of the molcTcules (in which case the rotational 
motion will not be effective), expressions for spectral 
densities will be obtained. The case in which the spin is 
not at the center of the spherical molecule does not alter 
the arguments of this section, and will be briefly discussed. 

The interaction hamiltonian can be expressed as 
(see Eq, III-24): 


«'»> ■ I z z 

3 q, P 


pCL qP T5,-(P+<l)> 

°P ^3 




(III-31) 


where the sum is over all j nuclei external to the molecule. 
The spectral densities Eq. (III-25) in this case will then 


become; 
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4-00 



PjP'' 


P'5" p-(P+4)(t) F-(P'+4')rt-x) 


(III-32) 


21ais expression is ?/ritten by neglecting the cross-correlation 
of dipolar interaction between different pairs* This assump- 
tion is made without any a priori justification [35] j for want 
of a better treatment. Purthei^ since the magnetic energy 
levels designated by different values are closely spaced 
compared to kT and since there is a very large number of ^ j 
nuclei to be summed over, it may be assumed that the neigh- 
bouring nuclei are evenly distributed amongst their ( 2S+1) 
levels [ref. 36, p. 233]. By doing so we are effectively 
treating the spin S as, part of the lattice'*' which is assumed 
to remain always in thermal equilibrium. This summation may 

then be carried over by replacing the product S? by its 

D <3 

expectation value 6^ ( |p }+l)S(S+l)/3. 

Using the correlation futiction for translational 
diffusion caloulated under "extreme narrowing"^ assumption 
[ref. 7, p. 300, and 35], the spectral densities (III-32) 
are obtained as: 



(III-33) 
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where N is the ntunher of spins (or molecules, for this case) 
per unit volume, ’^a'^ is the radius of the assumed spherical 
molecule and D is the transla,tion diffusion coefficient. 

The translation correlation time is then given hy = 
2a^/D. Equation (III-33) again represents a random field 
(compare Eq. III-23) having a mean square value; 

f = 2lE_ Yq (III-34) 

^ 5a^ lb 

Thus under ext reme-nanr owing approximation, and 

if the spin states of S can he assumed to he equally popula- 

z 

ted (high tempera,ture approximation), the intermole cular 
dipole-dipole interaction can*he approximated hy a random 
field. 

This conclusion remains unaltered if the spins are 
considered to he situated at a distance from the centers 
of the spherical molecules. In this case the average distance 
of closest approach between the spins becomes different from 
the above value (2a) and further both translational and 
rotational diffusion contribute to the time dependence of 
this interaction. These two are compensatory effects. The 
correlation function calculated for this model, for b ^ a/2, 
by Hubbard [35] leads to an expression similar to (III- 33) 
with a slightly modified value of mean square fields which 
are given by; 

f = Yj y| ‘^^Ci+0.233(^)^ + 0 . 15 (^)^ + .....] 

^ 5a^ 


(III-35) 
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Tlie second and third terms add to ‘less than 7 per cent of 
the first term for To = a/2» 


Correlations ; - If there ; U ■ more than one spin in the 
molecule, the interaction hamiltonian (III--31) will "become: 





(p+q.) 

(tjj 


(111-36) 


where ^i^ stands for all nuclei in the molecule. The various 
terms th^t will occur in the calculation of spectral 
densities will "be, in general, of the type; 



The terms for vdiich i = i'^'and j = j represent the 
completely correlated terns, included in the aloove calcula- 
tions.' The terms for which 3 ^ 3 ^ represent cross-terms 
"between dipolar interactions with external spins and may 
Toe assumed to Toe completely uncorrelated and talren as zero 
[35] (these terms are zero even when i = i'^ , but 3 7 ^ 3 ^ 
and dt is this part which v/as assumed zero in Eq_. (111-32)) » 
The terms for which i ^ i^ but 3 = 3 ''’ represent the field 
produced at the two spins i and i'' in the molecule by the 
same external spin 3 and represent the case of partial 
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correlation, lloggle [14] has performed a calculation of 
this term hy assuming only translational diffusion gnd 
neglecting the effects arising due to the spins heing not at 
the centers of assumed spherical molecules, Y/hich also leads 
to f(i) = f(i''). It is further assumed that the surface of 
closest approach is anellipsoid of revolution about the 
internuclear axis of the tvro spins i and i^ having a semi- 
major axis (r + h) and semiminor axis h, where 2r is the 
distance between the spins i and and h is the nominal 
distance of closest approach of spin j to any of the spins 
i (assumed equal for both) , The calculation showed that for 
h r, C = 1 and for r h, C -^0. 


(b) Scalar Coupling :- Using the hamiltonian (III-IO) and 
Eqs- (111-2,11) for the scalar coupling of a single spin ^ 
with a spin S. having fast time dependence in its spin 
states due to fast quadr^ipole relaxation, and using Eq. 
(III-21), the spectral densities for this case are obtained 


as: 


T-nv 


ioo 





-oo 




iCJy qT 

e 

( ( 1 1 ^ 1 + 1 ) 




<(^~^(t) S“^'^(t-T)^ 


av 


(III-58) 


Assuming Bloch equations for the time dependence of spin ^ 
(ref. 7, p. 310 and [3]), one obtains q = -q'' and : 
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jiiv ^ lql+1) 1 






(III-39) 


Using extreme narrov/ing (<jOqT^ ) <?C 1 where it is the 

correlation time for molecular reorientation viiicli causes 
the fluctuations in the electric field gradients at the site 
of spin S (this assumption has to he made for Bloch equations 
to hold for motion of spin S for S ") 1, [3]),vre obtain = 

/-o> / 4- 

Tp = T where t is the quadriipole relaxation of S. 

Equation (111-39) then represents a random field (compare 
Eq.III-22 and note (^j instead of in the denomi- 

nator) with a mean square value: 


4tc^J^S(S+1) 

5 


(III-40) 


Thus the only condition for scalar coupling to be 
approximated as a random field is '''' extreme narrov/ing f or 
molecular reorientation, which is generally true, _ It may, 
however, be noted that (CJ^ - is normally not 

less than unity and in fact quite often ( ^ 1 
is obtained. 

Correlations;- If there are number of I(i) spins coupled 
to one- S spin (having fast time dependence in its spin 

/-w 

states) then the random fields at various I(i) spins will 
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"be completely correlated (0 = 1 ) with, mean square fields 
'equal to 47i^ J? ^S( S+l)/3. 

lo 

In case there are more than one spin of the type 
8 ( 3 ) (having fast time dependence) the portions of the 
fields produced at the site of any spin I(i) hy different 
3 ( 3 ) spins may be assumed to be completely uncorrelated 
with each other. This then leads to a case of partial 
correlation between the fields at various l(i) spins with 
the extent of correlation governed by the coupling constants 


(c) Spin-Rotation Interaction ;- Using the hamiltonian 
(III-18) and Eq. (III-21) the spectral densities are 
obtained as: 


q-oo 



This expression has been evaluated by Hubbard [31], under 
the assumption that the characteristic time for the 
change of the components of angular momentum is very 
much different from '^2 molecular reorientation, as: 


..nv 






(2G^+G„)^2t-^ 
[l+(U}j qV]_)^] 




[l+(U)j 


( 111 - 42 ) 


3 
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T/here 


Ti2 t:3_ ^ ^2 


(III-43) 


and either <<( X 2 '^2 ^1* niimher 

of identical I spins on the same molecule and Icl is 
the Bolt2mann factor. In writing expression (III-42), 
no correlation hetv/een spins on different molecules and 
complete correlation "between spins on the same molecule 3Ye 
■ ass-umed [33]. 


Eg. (III-42) represents a sum of two random fields 
having different correlation times and different mean square 
fields. However if j "^12 ^ '’'I (III-42) gives 

a random field with mean square field as; 



(2C_j_+Cjj)^ + 2(Cjj 



= H 


Sf ( 2Cl+0p , 


3t' 


(III-44) 


IIIG TREATMEII OF SCAI.IR COUPIIHG AS A EELAXA.TIOH IIECHAIIISM 

?/hen a spin 5 nucleus is coupled by scalar coupling 
to a nucleus having spin 1, this coupling generally does 
not appear as a splitting due to the rapid quadrupole 
relaxation rate of the latter. Depending on the relaxation 
tiwie T of the quadrupole nucleus and the magnitude of this 

ml 

coupling a broadening in the spectrum of spin ^ nucleus is 
generally observed. 7i/hen this coupling is entirely ’^v/ashed 
out [37] (no residual splitting remains and only broadening 
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is observed) the scalar cotipling interaction is treated as a 
relaxation mechanism. In the following the conditions under 
which such a treatment of scalar coupling becomes valid are 
discussed. 

The manner in v\?hich the coupling of I (I = -|), with S 

./W 

(S ^ 1) affects the' single resonance spectrum of spin ^ can 
be treated by considering the total hamiltonian: 




+ 


(III-45) 


where 


ifll are zeeman terms of spins I and S 


respectively, '^£sc scalar coupling between them, 

is the interaction of the quadrupole moment of ^ 

with fluctuating electric field gradients and is responsible 

'IJ 


for relaxation of spin S, interaction of 

observing r.f, field with the spin system and 7^^(t) 
represents mechanisms of relaxation of I other than scalar 
coupling yath S. The equation of motion of the spin density- 
matrix o can then be written ass 


do 

dt 




Ca-Oo). 


( 111 - 46 ) 

This equation includes complete information regarding the 
dynamical behaviour of the spin system including spin S. In 
order to describe the relaxation effects on a transition 
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(m(i) -T^m(i'')) in the spectrim of spin I, one considers 
simultaneous existance of off-diagonal elements of 0 of the 
type |m(i), m^^ k — > jm(i^), lor all different values of 

m^, the eigenvalue, of S^, Purthermore^as may he seen from 
Eq_. (III-39) that ifi4n:^J^^^ (^j - ^3)^ ff the relaxa- 
tion time of spin S due to f® such that 


(bU - tyJo)^ yy Ij the terms like (l,S + I S , ) in I*S 

X D CJ, ' -j-. — <-« -j- ^ ^ 

will not contribute significantly to various relaxation 
matrix elements. In such a case only th& term Jl^^g 
retained in c- #ith such an approximation and setting 
= 0 in Eq. (IH-46) . an expression for 1=4- and 


S = 1, for the line shape of the spectrum of spin I, is 
obtained 
as [38]: 


obtained, with the relaxation time x of spin S as a parameter, 


I(x) 


2p 

2nJ 


(45+Tl®+5x®p^) 

225x^ ( 34x^ - 2x^ +4 ) +t]'^(x^-2x^+x^ ) 


j (III-47) 


vdiere p = lOmT^^J and x = A6vJ/2-n;J, in which is the 

difference between the laimor frequency of spin I and the - 
observing frequency. .Equation (III-47)» given by Abragam 
[7j p. 504], was also obtained through basically similar but 
slightly different arguments by Pople [39] following a 
general treatment of lineshape by Sack [40] and Anderson [41]. 
This equation describes the spectrum of spin 1 for the entire 
range of values of If <[<[ 1, Eq, (III-47) leads to: 


I(<^) = - 


2p 

2kJ 


45 

225x^+4p^ 


3/'8 2 _ -rZ 

XgJ 




) 


( ALO)2+(Stc^x 


(III-48) 
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a lorentzian having a half-width at half-intensity equal- to 

8%^-r^jV3. 


This is also the linewidth and line shape obtained 
when scalar coupling is treated purely as a mechanism of 
relaxation with the condition that ~ 

Prom the expression for spectral densities for scalar 
coupling mechanism Eq. (111-39), (assuming extreme narrowing 
for motion of spin S) one obtains for T-, and T, of spin I: 

J-iC ^ 


1__ 8Tt^J^S(S+l) 

^1 ^ i+(a?3- 


(III-49) 


and 


L 

Tr 




8(5+1) 


11 . 




+ x 


1 


(III-50) 


Prom Eq. (III-50) it can be seen that, under the condition 
^ 1. 8=1, the linewidth (l/T^) is 

the same as in (III -48) . 

Thus the conditions under which scalar coupling can 
be treated as a relaxation mechanism are that t]^ <^<C 1 and that 

V jL 

extreme-narrowing'*' holds for motion of spins 8. 

A-' 

Equation (111-49), shows that if the above mentioned 
conditions i.e, 4m^J 2 «(CJ^ - Co^)^ and (U)j - iA>g)"T^ y^ 1 
hold, the contribution of scalar coupling to l/T^ becomes 
negligible and in that case the scalar coupling mechanism 
serves only as a source for spin-spin relaxation. The contact 
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of tJie spin system with the lattice is therefore primarily 
through other interactions in the sample. 


HID GROSS-IEEIvIS BETWEEN VARIOUS INTERACTIONS 


In a sample various interactions mentioned in the 
previous sections often simultaneously co-exist and it is 
only their relative importance that differs from sample to 
sample. In such a case there can exist cross-terms Between 
different interactions in the calculation of spectral 
densities such as: 



(III-51) 


Besides cross-terms within an interaction. The cross-teims 
within an interaction have already Been discussed under the 
title '^correlations'^ in each case and in many cases found 
significantly important. The cross-terms Between different 
interactions normally vanish Because of either very different 
orders of magnitudes of time dependences involved in their 
hamiltonians,iu which case the two parts can Be separately 
averaged giving zero for the averages, or difference in 
transformation laws for each of the hamiltonians yielding 
zero on averaging [ref. p. 509], Specifically the • 

cross-terms Between scalar coupling and either of the dipole- 
dipole interactions or the random fields vanishes Because of 
different orders of time-dependences involved. 
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17 - EXPESIMEUJOIAI 

I 7 A FEEQUEHCY SWEEP DOUBLE EESOITAUGE SET-UP 

!• Introduction 

In homonuclear douUle resonance experiments tlie 
irradiated and' the observed transitions Delong to the same 
nuclear species, and therefore the frequencies of the 
irradiating and the ohserving r.f , fields differ hy an audio- 
frequency. In such experiments the second r.f, field is 
produced hy modulating the magnetic field with the appropriate 
audio frequency. Douhle resonance experiments can, in 
principle, he performed both hy changing the magnetic field 
and keeping the frequency fixed (field sweep operation) or 
hjj' keeping the magnetic field fixed and changing the frequency 
(frequency sv;eep operation). The frequency sweep experiments 
are convenient for relaxation studies and are used in the 
present vrork. 

The conventional high resolution n.m.r. spectro- 
meters, like the Yarian IIP-lOO, use field sweep operation. 

This instrument has been modified to perform frequency sweep 
single and douhle resonance experiments, hy using field- 
frequency-locking control loop [25,25]. Single resonance 
spectra have been obtained without loss of resolution and 
with improved stability, compared to field sweep operation. 

The modification, however, does not disturb the field sweep 
operation. The details of this modification are described in 
this section. 
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2» Basic Principla 

111 tlie first stage of this experiment the magnetic 
field is held fixed at a particular value. Ihe experimental 
sample is added^j/tha small anomit of tetramethylsiiane 
(SliS) v/hich serves "as a source of the locking signal. The 
magnetic field is modulated hy an audio freq,uency f^^, Pig. 
(17-1), and is brought ill resonance \?ith the first lorr-field 
sideband of TIIS. The output of the r.f, receiver is then 
phase detected , through a loclr-in-ampllf.ie.r. which is 
referenced by the same audio oscillator and has amplifiers 
tun^sd to frequency f^. The phase difference between the 
signal and the reference of the lock-in-amplifier is so 
adjusted that, as one sweeps through the resonance, the d.c. 
output of the lock-in-amplifier is obtained in the shape of an 
error signal. This error signal is fed to the flux-stabiliiser 
such that it gives an appropriate correction voltage whenever 
the field drifts in either direction from its locked value. 

The output of the r.f. receiver is also monitored on a scope, 
after amplification through an amplifier tuned- to frequency 
f^. A continuous steady display of a signal at f^ indicates 
a locked field. 

The single resonance spectrum is obtained by modu- 
lating the magnetic field by a second audio oscilla,tor and 
STiTeeping its frequency (f 2 ) by a synchronous motor, ITlien 
^2 ~ ^ 6, where 6 is the frequency difference between the 



Mg, IV~1. Block diagram of the experimental arrangement 
for a frequency sweep homonucloar double 
resonance set-up. 
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signal to le observed and IMS, tiie sideband of the signal mil 
come in resonance with the loclied field. The output of the 
receiver which now contains a component at fg also, is fed to 
another lock-in- amplifier which is referenced by the oscilla- 
t car gene rating fg. This lock-in- amplifier is a broad-band 
amplifier. The -output of the lock-in- amplifier is fed to a 
recorder. The sweep oscillator has a linea,r frequency dial 
so that the different parts of the spectrum are all scanned 
at the same rate and that the spectra are obtained with a 
linear frequency scale. 

The irradiation of the spectrum is achieved by 
modtilating the magnetic field with a third a,udio- oscillator 
(f^) such that f^ = f^ - 6^^, where 6 '^ is the frequency 
difference between the line to be irradiated and TMS. The 
amplitude of this oscillator is controlled by a calibrated 
attenuator. 

5* Details of the Experim ental A rrangemeni; 

The high resolution n.m.r, spectrometer HR-100 

consists of a 12 electromagnet (max, field 23.5 x 10^ gauss) 

havir^ a current regulated power supply. A flux stabilizer 

V3506 controls the time variations in the magnetic field while 

a field homogeneity control unit smoothers out any spatial 

variations of the field in the region of the sample. Crystal 

controlled r.f. transmitter and receiver V4311 operating at 

100 MHz are used for proton resonance. The homogeneity and 

8 

the stability of the magnetic field are^l in 10 . The 
resolution of the instrimaent is '—'0.5 Hz. A linear sweep 



80 


I 

unit Y4352 giving a saT'-^tootli is used for scope xDresentation 
of the signal and a slovr sv;eep unit V3507 which operates 
in conjunction with the flue-: stabilizer for slo\7 recording 
of the spectrun in field s^veep mode. 

(a) lock Channel ;- Ihe purpose of this channel is to hold 
the magnetic field fixed at a particular value, She 
2 IvEz audio-oscillcutor and lock- in-amplifier available in 
Varian V3521A integrator, have been used for this purpose 
Pig. (IV-2) . The output of the audio-oscillator is attenuated 
by a GR 1450-TB attenuator, which serves as an isolation 
stage, before feeding to the sv/eep coils. Since the lock 
signal is derived from side-band resonance of TIIS, the 
integrator is operated in side-band mode. The receiver 
output at j314 . on Y4311> v;hich nov/ contains an audio compo- 
nent at 2 ICHz, is fed to the signal channel (referenced 
internally) through an attenuator. The phase difference 
between the reference and the signal can be adjusted with a 
front panel control. The output of the lock-in-amplifier, 
which is the error signal^ is taken from the ^^recorder output 
of the integrator and fed to the flux stabilizer through a 
switch. This error signal may be fed to pin 5 of J1 on the 
sloTf sweep unit. The amplitude and phase of the error signal 
are adjusted for maximum stability of the lock. The stability 
is monitored by feeding the output of the receiver at J314 
to a scope, after amplification through a GR 1232A tuned 
amplifier tuned to f^. An amplitude of the order of 100 
millivolts is usually needed for the error signal. 
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Fig. IV-2. Blook diagram of the look-ohajmel. Iho 

oomponents shovm mthin the hrokcn lines are 
available in the Varian integrator V3521. 
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(b) Observation. Chaimel s-- A. OR 1304-B, beat-freqnency 
audio-oscillator, equipped with, a linear increment dial 
(range 100 Hz) is used in tliis channel, fig. (IY-3) . 'Ihe 
output of this oscillator is fed to the s?/eep coils at the 
modulation input jack on linear sweep unit through a HP 350-D 
attenuator. She lock-in-amplifier used in this channel is an 
" EMC ''(Electronics Missiles and Communications Inc. Mount 
Yemon, IJSA) model EJB. Since the imput impedance of the 
reference channel of this lock- in-amplifier is 100 Eil,, the 
reference signal from beat-frequency oscillator is fed 
through a imp e dene e matching net work. 

Ihe signal channel of the MC lock-in-amplifier 
contains a negative feed back pentode amplifier which is 
heavily degenerative except at null frequency of a t\^/in-tee 
net work, placed in the negative feed back circuit. Since 
the frequency is swept in the present experiment, this lock--, 
in-amplifier is to be operated over a range of frequencies, 
Ihese twin-tee •’’s were, therefore, eliminated from the circuit 
and the screen grid of the pentode amplifier grounded, and 
as there is no degenerative feed back at any frequency the 
circuit acts as a broad band audio-amplifier. 

Ihe output of this lock-in-amplifier is obtained 
through a low frequency band pass filter and fed directly 
to a G-14 recorder. The phases of the frequency sweep 
spectra are adjusted by adjusting the reference phase of the 
EMC lock- in-amplifier. 
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Fig 


I?- 3 , Block diagram of fhe observation channel. 

/in EMC, model EJBj, look-in-amplifier and a 
GR 1304-B beat-freq[uency audio oscillator 


are used in this chjiinnel. 




OBSERVATIOK’ CHANNEL 
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(c) Irradiation Qliaimel ;- A HP 205-AG- audio-oscillator lias 
been used for irradiation. Iliis oscillator has a maximum 
output povrer of 5 watts and a calibrated attenuator in its 
output, The output of this oscillator is fed at the 
modulation input jach on linear sweep unit. 

A HP 521C frequency counter is used for measuring 
frequencies of various audio oscillators. 

(<i) Pield-frequency Control Panel: - A control panel 
essentially consisting of various switches, for a quick change 
from field-sweep operation to frequency sweep operation and 
viceversa has been installed just below the oscilloscope. 

One switch. Pig. (IY-4), changes the input to the oscilloscope 
Y between integrator ''scope output and the tuned amplifier 
( 2 KHz) output for monitoring the lock. Another switch 
changes the recorder from integrator ’^recorder-output* 
to ^^recodrder-output'^^ , One oii/off . switch 

is used for the error signal to the flux stabilizer. The 
amplitude of the error signal can be adjusted by the poten- 
tiometer Tune lock 

If all the switches in this panel are flipped 
down the spectrometer is in the field-sweep mode. 

(e) Operational Details ;- The magnetic field is adjusted 
for best homogeneity and resolution in the field sweep mode 
and the low field side band of TMS ohtelned in the shape of 
an error signal in this mode. At the mid point of this 
signal the sweep is taken off and the error signal fed to the 
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Fig* IY-4. Diagrsjn sliomng the functions of various 
switches on the field-freq_'uency control 
panel , 




field-frequency CONTROL UNIT 
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flux stabilizer, _ Ibe sign of the error signal is chosen by 
trial and error by noting the locking or antilocking of 
the field. Ihe lock is made stable by adjusting the ampli- 
tude of modulation and the amplitude of the error signal 
being fed to the flux stabilizer. A slight adjustment of 
the r.f, reference phase at V4311j has also been found 
useful for stabilizing the lock. The stability of the lock 
can be tested by deliberately trying to drift the field 
slovvTly and noting that the error signal compensates for it. 
Small variations in homogeneity are adjusted by maximising 
the output of f^ at J314, v/hich is monitored on the scope, 

Eor the purpose of scanning the spectra the linear 
increment dial of the beat-frequency oscillator has been 
coupled to a synchronous motor through gears. Different 
driving speeds are obtained by changing the gear ratios. 

The amplitude of the beat-frequency oscillator is 
usually kept large for the following reasons. As the input 

impedence of the reference channel of the EMC lock-in- 

the 

amplifier is 100 the output of A.beat-frequency oscillator 

is first dropped across a 1 K-^ resistance to match the 
output impedence of the oscillator and then dropped across 
a 100 Kru resistance to match the input impedence of the 
reference channel of EMO (Eig, IV-3) , The two resistances 
are isolated by a 100 KjTl potentiometer. The modulation 
amplitude of f 2 and the signal to the reference are adjusted 
respectively by the HP 350-D attenuator and the 100 K iX 



90 


potentiometer, A large amplitude of the 'beat-freq.uenoy 

oscillator then ensures a high attenuator setting and a 
‘tbs- 

high value ouaIOO KjIl. potentiometer, A high value ofAlOO KjTL 
potentiometer avoids any mismatching of the reference input 
impedence. ' 

Ihe frequency of irradiation is made to exactly 
coincide with a single resonance transition in -the following 
manner. In the absence of irradiation the observing 
frequency f2, is swept to the peak of the line to be irra- 
diated and the motor is then switched off. The amplitude of 
irradiation is then made fairly low (of the order of the 
amplitude of the observing f I eld ) and its frequency 

adjusted as near to f2 as is possible with a frequency 
counter (within fraction: of a Hz), At this stage the output 
of the MG lock- in-amplifier will beat at the difference 
frequency of these two oscillators. This beat can be 
monitored either at output meter of the EMC lock-in- amplifier 
or on the recorder. The frequency of irradiation is now 
carefully adjusted by a fine control to reduce the beat 
frequency to zero, at which the EMC lock- in-amplifier will 
show a large d.c, amplitude as indicated by a large deflec- 
tion of the output meter. 

The irradiation oscillator is normally operated at 
full power and various amplitudes of irradiationAobtained 
by using the calibrated attenuator. This facilitates (if 
all the other parameter are kept constant) a calibration of 
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complete range of strengths of irradiation once one of the 
values is calibrated. 

Performance 

A comparison between the field-'sweep single 
resonance proton spectrum of 2-bromothiaz:ole and its freq- 
uency-sv/eep spectrum obtained with the above set-up is shown 
in Pig. (IY-5) . The resolution in the frequency sweep 
operation is of the same order as the field sweep operation, 

IVB SAf.'IPlE PREPARATION 

Commercial samples have been used in this work. 

The samples were purified, by recry stalizing in case of 
solid samples and redistillation in case of liquid samples, 
before sealing in 5 m.m. outer diameter n.m.r. sample tubes, 
A few drops of tetramethylsilane was added in each sample 
to provide the locking signal. The samples were degassed 
under vacuum by a series of freeze -pump -thaw cycles [14,22] 
and finally sealed in vacuum. 
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5’ig. IY-5. The single resonance proton spectrum of 
2-1) romothi azole (a) in the field swoop 
operation and (h) in the frequency sv/eep 
operation obtained with the present set-up. 
The broader doublet is assigned to proton 
(position 4). The chemi’cal shift between 
the tv/o groups is | *1^ - 'Pp | ” 18.8 + 0,1 
Kz and the coupling constant = 3.55 + 

0.05 Hz. 
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Chapter V 

DOUBLE BESOMCE SgUI)Y OF PROTOh EBLAXAglOK II A 
SYIffiETRIOAL THESE SPIN SYSTEM (ABg) 
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V - DOUBLE BESOHMGB STuLI OP- PRQgOIT BELAXASIOIT IH A 
' SYBMETSICii THESE SPIN SYSTEIa (ABg) 

VA INTRODUCTION 

In -fehis chapter the results of a proton-proton 
double resonance study of spin relaxation in a symmetrical 
three spin system (AB 2 ) formed by the ring protons of 2,6- 
dibromoaniline are presented. The sample consisted of a 
20 per cent by weight solution of 2,6-dibromoaniline in 
COl^. Double resonance spectra exhibit significant relaxa- 
tion effects over a wide range of strengths of irradiation 
(from 0.1 Hz to 9.0 Hz). The effect of magnetic field 
inhomogeneity on the linewidths of double resonance transi- 
tions has been explicitly considered in this analysis and " 
is found to be of considerable importance in deriving infor- 
mation on the relaxation processes. Tvro relaxation mecha- 
nisms s(i) external isotropic random fields and (ii) internal 
dipole-dipole interactions between the ring protons along 
with possible correlations between them, have been considered 

If 

for detailed calculations. The calculations for high 

\\ u 

strengths and low strengths of irradiation have been 

performed by the methods outlined in Chapter II. A detailed 

5)3ectra 

comparison of the calculated’ A for both' high and '’low’ 
strengths of irradiation with experiment is found 

to be essential for distinguishing between the contributions 
of the different mechanisms in this study. 
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VB SINGLE RESONAITOE SPECTRUlv! 

PjS-Dibromoaniline Belongs to a symmetry point 
group 02 ^. and its ring protons form an type of spin 
system (in tlie notation of Pople et.al. [2]). Tlie resonance 
of proton"!*^ (in the para position with respect to NH 2 
group) is expected to appear at a higher field than that of 
protons ' B*" (meta to NH 2 group) This gives 1)^ where 

1 ) is resonance frequency. The Basic-functions, diagonal and 
off-diagonal matrix elements of the spin hamiltonian and the 
transition frequencies and intensities for an ABg spin system 
are given in TaBle(V-l). The diagonal and off-diagonal matrix 
elements of double resonance hamiltonian in the rotating 
coordinate are also indicated in this table. In single 
resonance there are six symmetric and two antisymmetric spin 
functions, whose frequencies depend on two parameters viz. 

J^B» spin-spin coupling constant, and - Vg “the 

chemical shift Between the ^ a " and ' protons. 

The single resonance spectrum of the molecule is 
shown in (V-1) ♦ The spectrum has four well resolved A 

transitions and two closely' spaced doublets in B* transitions. 
The small separation in the two '''b^ doublets has been 
determined in the field sweep operation by wiggle-heat 

- i " mm . — — ■ 

NH 2 is an electron donating group and produces greater 
chemical shielding at the site of the para proton than that . 
of the meta proton in monosubstitute d benzenes, while 
bromine is inactive in this respect [ref. 2, p. 259]. 



Basic product functions, diagonal and off-diagonal matrix elements, wave -functions 
and energy values in single resonanoet of an MS>2 spin system for which jVoB}> [Vq^^I . 
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Pig. V-1. Prequency sweep single resonance spoctrura 
of tiie ring protons of 2,6'-dil)romo<aniline , 
which form an ABg typo of spin system. 'L'ho 
chemical shift hetween the two groups is 
I Ml " . = 89.95 Hs and the coupling 

constant is = 7.95 Hs. 






-■? .'/‘■■''rV?A' - 
Wllfewic* 
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method [ref, 2, p. 76], In this method if the separation 
between the lines in a doublet is J in cycles per sec., the 
time between successive beat maxima is equal to 1/J sec. 

The values obtained for this case are 0.3 + 0.1 Hz and 

0. 4 + 0.1 Hz respectively for high and low, field doublets. 

The analysis of the single resonance spectrum is straight 
forward and the parameters obtained are j I = 

89.95 + 0.1 Hz and = 7.95 + 0.1 Hz, 

The different transitions are shown schematically 
with the help of an energy level diagram in Hig. (V-2) 
and are nvimbered in the order of increasing frequency. Lines 
1 to 4 are A transitions and lines 5 to 8 are B transi- 
tions. line 3 is the transition between the antisymmetric 
states 7 — ^ 8. 

VC DOUBLE EESOlIiUTCE PECTRA 

1. General Eeatures 

Double resonance spectra were obtained by irra- 
diating each one of the lines of the single resonance 
spectrim, for values of the strength of .irradiation, 
varying between 0.1 Hz and 9.0 Hz. Some of the typical 
double resonance spectra are shown in Pigs. (V-3 to 7). 

The spectra with'*low'"’^strengths of irradiation do not show 
any splittings but do show marked intensity changes compared 
to single resonance, while the spectra with high strengths 
of irradiation show all the features of double resonance 
spectra viz, line splittings, line shifts, intensity and 
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Pig. 1 - 2 . Energy level diagram of fxn. AlDp spin nyntoin 

for which '"2^ ^ H* spin functions nol/C'd 

against different energy levels arc tluj 
respective spin product functions. '.I’h>; 
various transitions are £ilso sliovm. sciiomatl- 
cally with heights proportional to tlio rjiugle 
resonance intensities, 'fho linos 1 to 4 aiv 
'l'** transitions and the linv.'s 5 to 8 arc 
‘b"" transitions. Line 3 is the traasi tiion 
between antisymmetric states 7 and 8. 
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linewidth clianges. The values of Vg for spectra which 
show splittings were determined in the following manner. 

The hamiltonian in the rotating coordinate (E<1« 

is diagonalized by a computer for various values of .V 2 f.or 
any particular line of irradiation and the calculated line 
positions in the double resonance spectra plotted against 

These plots are then compared with the observed spectra, 
for irradiation of that line [11]. This leads to a deter- 
mination of V 2 to an accuracy of/-^+ 0.05 Hz, For cases of 
weak irradiation, low values of Vg were produced in the 
experiment by using a calibrated attenuator, and therefore 
for these spectra the values of Vg were calculated by 
multiplying the values determined for strong irradiation 

cases, by the respective attenuation factors, 

\\ > 

The spectra for low amplitudes of irradiation 
exhibit the following general features. Figure (V-3( a)) 
shows the spectrum of A group while line 6 of , B group 
is irradiated by an amplitude of 0,1 Hz, This produces 
a decrease in the intensity of line 1 and an increase in 
that of line 2, The corresponding spectrum for irradiation 
of line 7 (Fig. V-3(b)) sho?/s a reverse effect in the 
intensities of lines 1 and 2, Both these features can be 
qualitatively understood by referring to the energy level 
diagram (Fig. Y-2) and noting that the irradiation of a 
pair of levels causes an increase in the population of 
higher level and a depletion in that of the lower level. 
Therefore for irradiation of line 6 (6 4 transition) the 
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Fig . V- 


OlD served and calculated double rosonanco 

\ sf 

spectra of A group in 2, 6-dibromomiilinc' , 
while (a) lino 6 and (b) line 7 of B group 
irradiated with an amplitude of 0,1 llz, i'hi; 
calculated spectra are for random field 
mechanism (mechanism (i)) with f(A)/;r(B) = 


rad. /sec. 


= 0.0 and 2f(B)T:^ = 0.2 

G 
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intensity of line 1 (6 5 transition) should decrease 

while that of line 2 (4 3 transition) should increase, 

in conformity with the observation (Pig. V-3(a))« These 

the 

changes are reflected in^ intensities of all those transi- 
tions in the spectrum which have an energy level in common 
with the pair being irradiated"^. Similar intensity changes 
have also been observed for weah irradiation of other lines 
except line 3. Typical cases of "b*' spectra when lines in 
A group are irradiated wl±h - 0.1 Hz are shown in 
Pig. (V-4). 

A weak irradiation of line 3, which is the transi- 
tion between the antisymmetric states? does not produce 
intensity changes in other lines'i*"^ and furthermore in any 
of the observed double resonance spectra (Pigs. V-3,5 and 6) 


f Intensity changes, smaller than the above, can also occur 
in the other transitions which are connected by relaxation 
processes to the irradiated pair. These changes are not. 
obvious from an inspection of the spectrum but may become 
significant for understanding the spectrum' on a quantita- 
tive . basis. 

■|-|- Intensity changes and splittings occur for irradiation 
of this line for large values of V 2 , due to the fact 
that the symmetric states belong to the neighbouring 
transitions (particularly line 2 (4 3 transition) , 

get perturbed by the irradiation. 
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Pig. V- 


Olservod and calcnLatod doullc resonraine 
spectra of ‘b* group in 2 , 6-dil)romof:iniiinu, 
.vhile (a) line 1 .md (b) line 4 oT 'a' group 
are irradiated with an amplitude of 0.1 Hb. 
Iho caloulated spectra arc for random tiuld 

mechanism V7ith f(A)/f(B) = 2.0, 

G ^ = 0.0 ,and 2f(B)'rQ = 0.2 rad./seo. 
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Observed and calculated double rosoriatice 
spectra of *a' group in 2 , e-dltromorinlllno , 
while (a) line 6 aiid (h) lino 7 of ‘b' group 
are ixx’adiated witli an, ajnplitude o.fc 2.«- Hz. 
Ihc calculated spectra aje for random field 
meciiauism vath f(A)/f(B) = 2,0, ~ 

and = 0.0. 
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Fig. V-6. Observed and oalciilateddoiiblo resonance 

spectra of'^A* gronp in 2, 6~dibromo aniline, 

« I 

while (a) line 6 and (b) line 7 of B group 
are irradiated with an ajnplitudo of 6,7'J IIz. 
The calculated spectra arc for random field 
mechanism with f(A)/f(B) = 2.0, 0^^ == 1.0 

and = 0.0. 


A 
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Fig, y-7. Observed and calculated double rono'njsiac; 

spectra of *B* group in 2, 6~dibroinoanil:hif,; , 
while (a) line 1 .and '(b) line 4 of A." group 
are irradiated by amplutiidos of 5.25 Hz and 
9.25 Hz respectively, i’he calculated speeln'.a 
are for random field mechanism with f{A)/r(B) 
= 2.0, =1.0 ,and = 0.0. 
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tlie intensity of this line is not changed'^. Ihis indicates _ 
that the relaxation hamiltonian does not connect the symmetric 
and antisymmetric states. 


2. Linewidths 


The double resonance spectra also exhibit significant 
linewidth variations, For example the spectrum of 'a* proton 

I s 

obtained when line 6 of the spectrum of B proton is irra- 
diated by an amplitude of 2.2 Hz (Fig. ■V-5(a)), has tvro trip- 
lets nearly symmetrical in frequency with respect to center, 
but having considerable difference in their linewidths. This 
asymmetry gets exactly reversed when line 7 is irradiated by 
the same amplitude (2.2 Hz) (Fig. 7-5(b)), It was pointed 
out in Sec. (IIB-2) that the inhomogeneity in the static 
magnetic field can contribute differently to various transi- 
tions in the double resonance spectrum and that this contri- 
bution to any transition a — s may easily be obtained 

from the plot of versus for weakly coupled spin 

systems. For making this estimate the present system may be 
considered as a weakly coupled system to a good approximation 
(J/ “■ Vjl 0.09). A plot of the frequencies of the k' 

transitions versus A^ = ( /2% for 

= 2.2 Hz, is shown in Fig. (7-8). These plots are 

P R 

constructed by diagonalizing -Mq with the help of the 


-j- line 3/ therefore serves the useful purpose of being an 
intensity reference in the double resonance spectra. 
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Pig* V-8, 


Plots of oolcolatod doublo roso^oT.oo :|■.•o.,uoo- 
oios of ’a' tranal-tions versus truqueiie.v 

„ 4..:^^ of B niiolox. 

offset of irraaiation o.i. • 

noted for the AB,, Bpin GyBtem 

plots are o-aLoalatoa i- 

^ = 2.2 Hs. Tlio nmidorlni-, of the s-lui .e 

oorresponding to tire transitions, i-' •««■> 

that for largo positive values of Ap and 
= 0.0, those states Irooomc the slnfilo ^ 
resonance states as indicated in Kg. 
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computer and. de'ijemjx-ing the corresponding transition freg_uen- 
cies for different values of A^, keeping Vg constant. In 
joining the calculated points to construct these plots 
certain non-crossing rules [42] should he oheyedt The tran- 
sition frequencies in a frequency- sv/eep spectrum for any 
value of A-g are given hy the points of intersection of an 
ordinate^ drawn passing through that point, with the contours 
of lig. (V-8). . According to Bq, (11-63) the change in the 
contribution of the inliomogeneity to the width of a particular 
transition is then determined hy the slope of the corres- 
ponding contour at the point of intersection. It may he 
seen from Big. (Y-8) that the frequency of the transition 
between the antisymmetric states 7 8 is not changed^^due 

to a change in Ag, giving a zero slope for this line. This 
is a consequence of the fact that the irradiation does not 
mix states of different symmetry. This line (line 3) will 
then have the same contribution of inhomogeneity to the line- 
width as in single resonance. The frequency of the transition 
3 4, v/hich is weak in intensity, changes very little with 

A-q, since the coupling between A and B nuclei is not very 
strong in this case, ■ In a wealcly coupled spin system ^ 2 » 
this line merges with 7 — » 8 [43] and does not show frequency 


An operational statement of the non-crossing rule is; two 
contours belonging totwo transitions which have an energy 
level in common Virill not cross in the presence of a 
p e rtunb a ti on , 
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shifts. However^ in a strongly coupled spin system it 
should sho?; frequency shifts. 

Irradiation of line 6 corresponds to Ag = 3.6 Hz. 
Hrom the slopes of different contours in Hig. (V-8)^ at this 
value of Ag^it can he seen that for transitions 5 — 2, 

5—^4 and 5-^6, forming the low-fre quency triplet, the 
inhomogeneity broadening is mostly reduced while it is 
increased for transitions 1 — > 2, 1 —^ 4 and 1—^6, forming 
the high-frequency triplet in the same spectrum, in agreement 
with observation (see Big. ¥-5 (a)). (Through a similar argu- 
ment it is readily seen that this asymmetry is reversed when 
line 7 is irradiated, which corresponds to Ag = -3.9 Hz 
(see Hig. 5(b)). It may further be seen from this plot that 
at these ordinates the slope is nearly, +1 for various tran- 
sitions for which there is a increase and -1 for those for 
which there is an recovery in the inhomogeneity contributions. 
This gives an inhomogeneity contribution of 2 yAHq to the 
linewidth of the former and nearly zero for the latter. Hor 
transitions 3—^6 and 3—^2 the slopes at these ordinates 
are almost zero (Hig. V-8) giving the contribution of inhomo- 
geneity to the linevddths of these transitions same as in the 
single resonance. Mahing use of these observations, and 
the observed line widths for irradiation of lines 6 and 7 with 
Vg =2.2 Hz, and the arguments of Sec. IIB-2(c), a value of 
(yAHq)/tc is obtained to be about between 0.15 Hz - 0.18 Hz 
(full- width at half height). 



121 


3 # Galoulation of Eelazation B ffeots 

The theoretical double resonance spectra were calcu- 
lated according to the theory sketched in Chapter II for two 
relaxation mechanisms, (i ) isotropic external random fields 
and (ii) internal dipole-dipole interactions, the interaction 
hamiltonians and expressions for correlation constants for 
which have been given in Sec, IIIA, The spectra for large 
values of the strength of irradiation ( fYHgj z/’ \aa'^ ^ 

calculated using the Bloch approximation. Bor'lowXtrengths 
of irradiation, in which only intensity changes were observed 
and no line splittings and shifts occur, a solution for OC 
is obtained in the laboratory frame. 

In the calculations for mechanism (i) the mean- 
square values of the fields at the different nuclei, the 
correlation between these fields and the value of the 

correlation time^ are the unknovni factors in the absolute 
values of the elements of the relaxation matrix, For 
mechanism (ii) . x is the only -unknovra. fantor. 

G 

For mechanism (i), a fit with the observed spectra 
for strong irradiation can be obtained by performing the 
calculations for different valties of f(i)/f(i'^) (see Eq. 

III-4) and (see Eq. III-5), Tt has been pointed out in 

the previous section that the behavior of line 3, in the 
double resonance spectra obtained under different conditions, 
indicates that the relaxation mechanism In the mole exile does 
not connect symmetric and antisymmetric states. The 
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relaxation mechanism, tiierefore? seems to preserve the eq_ni- 

valence of the two b' nuclei implied by the spin hamiltonian. 

For rel.axation due 'to external random fields, the equivalence 

of the nuclei will be maintained if the random fields at their 

sites are such that (t) = (t). This leads to the 

■“l ^2 

condition that the random fields at B^ and Bg are equal and 
completely correlated [16], On the basis of this observation 
the calculations for mechanism (i) were performed by setting 
f(B^) = fCBg) = f(B) and 0^ ^ =1, and varying only the tv/o 

quantities f(A)/f(B) and C^^for high- field spectra"^ - 

Por mechanism (ii) the factors •g■(3cos® ) which 
determine the correlations betvreen the different pairs of 
spins (see Eq. III-9) are calculated assuming the geometric^ 
arrangement of the three spins in 2 , 6 -dibromoaniline to be 
the same as in an unsubstituted benzene molecule. The angle 
■"'tjjjj'' in benzene is 30 *^, ■ 30 ° and 120 °, giving 0.625, 0,625 
and -0.125 for correlation between (iB^, B^B^), (AB 2 , 
and 182 ) respectively. The distance between the protons 

is obtained (from the geometry of benzene) as 2.48 ^between 


+ It has been noticed in the computation that if f(B^) is 
set equal to f(B 2 ), a variation in not produce 

appreciable deviation in the intensity of line 3 from, its 
single resonance intensity. The relative intensity and the 
linewidths of other transitions in the double resonance 
spectruim axe also insensitive to this constant as will be 
shown by a later table (see Table V-5)« Thus even though 
Cb]_b 2 equal to unity to confiim to the theoretical 

requirements on the complete equivalence of the ‘‘B' nuclei, 
the experimental results do not allow an accurate estimate 
of this constant. 
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A aJid B protons and 4.29 A* Between B protons. 

For every set of parameters desoriBing tiie relaxa- 
tion process the ^ matrix is oBtained and suBstituted in 
Eq_. (11-59) to oBtain the peaJc signal intensity. Ail the 
other quantities in th.e numerator of tho,t expression are 
oBtained By diagonalizing ^ and oalcul opting the matrix 

v- O 

elements of I_^(i) and cr^ in that Basis. 'It should Be noted 
that the denominator in this expression contains the total 
linev/idth including contriButions from inhomogeneity and 
instability as discussed eahlier, This amounts to substi- 
tuting for i'fK.UJ) ^ in Eq. (11-59) j'^st the value of the 

in 

observed linewidth to oBtain tlie value of S ? 


VD DISCUSSION 


The caJ-oulated relative intensities in the double 
resonance spectrum for mechanism (i) with three different 
values 1.5,2#Q and 2.5 for f(A)'/f (B) and -G^=0.0 for all the 
three cases, and those for mechanism (ii) are given in 
Table (V- 2) along mth the observed relative inensities for 
the cases of irradia-tion of line 6 and 7 mth Vg = 2.2 Hz 
(see Fig. V-5) and for irradiation of line 1 with Vg = 5.25 
Hz (see Fig. V-7(a)). The value of given above is the 
-final value chosen. Postponing the discussion on C^ for 
the present, a. comparison of the theoretical intensities 
for mechanism (i) with the observed values shows that 
f(a)/f(B) = 2.0 gives the Best fit with the experiment. 
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Table V-2 

Obse rve d. ajid caloiilated relative intensities for irradiation of 
lines 6 and 7 ^'dtii an ainplitude of 2,2 Hz and . of line 1 ?ri.th an 
amplitude of 5.25 Hz. Ihe random field calculations are for 
C 31 B 2 - 1*0 ctnl = 0.0. For irradiation of lines 6 and 7, 
tiiG intensities a,re expressed relative to tlie intensity of 
line 3 (S.H. 5) -as unity. 


Irradia- 
tion at 

S.lt 

Observed 

relative 

intensity 

It 

Calculated relative intensity ^ 

Random field 

Dipole- 

dipole 



1.5 

f(A)/f(B) 

2.0 

2.5, 

1 

2 

3 


5 

6 



1 

0.33 

0.31 

0.30 

0.29 

0.28 ■ 


2 

0.63 

0.67 

0.64 

0.62 

0.65 


3 

0.42 

0.44 

0.42 

0.40 

0.40 


4 

0.54 

0.53 

0.53 

0.53 

0.53 

Line 6 


5 

1.00 

1.00 

1.00 

1.00 

1.00 

Vo=2.2 Hz 

C . 

6 

- 

0.04 

0.04 

0.04 

0.06 


7 

0.46 

0.44 

0.44 

0.44 

0.45 


8 

0.20 

0.17 

0.20 

0.22 

0.19 


9 

0.28 

0.26 

0 

* 

0 

0.33 

0.29 


10 

0.09 

0.06 

0.07 

0.08 

0.07 


1 

0.23 

0.11 

0.12 

0.14 

0.11 


2 

0.37 

0.29 

0.34 

0.37 

0.36 

L ine 7 

3 

0.22 

0.17 

0.20 

0.22 

0.15 

V2=2.2 Hz 

4 

0.55 

0.47 

0.46 

0.46 

0.51 


5 1.00 1.00 1.00 1.00 1.00 


0.04 0.04 0.04 0.04 


6 
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Table V-2(ooiitd-) 


*1 

JU 

.2' 

3 " 

4 

5 

5“ 

7 



7 

0.55 

0.49 

0.50 

6.50 

0.51 



8 

0.33 

0.41 

0.39 

0.37 

0.41 



9 

0.50 

0.47 

0.45 

0.43 

0.50 



10 

0.20 

0.14 

0.14 

0.13 

0.14 



1 

- 

0.02 

0.02 

0.02 

0.02 



2 

0.26 

0.26 

0.27 

0.27 

0.28 



3 

1.03 

1.03 

1.03 

1.03 

0.92 


Line 1 

4 

0.88 

0.95 

0.97 

0.98 

1.20 


V2=5.25Hz 

5 

0.78 

0.66 

0.63 

0.62 

0.65 



6 

1.00 

1.00 

1.00 

1.00 

1.00 



7 

0.25 

0.16 

0.16 

0,16 

0.23 



8 

- 

0.01 

0.01, 

0.01 

0.00 



f Serial rfumber of the double resonance lines in a decreasing 
order of frequency (left to right in the speotrum). 

■^■^These values are averages over several spectra (5-8) and are 
accurate approximately within 3 per cent of the values given. 

Calculated relative intensities have been obtained by dividing 
the calculated integrated intensities by observed linewidths. 





126 


Ihe oalcxiLated intensities for f(A)/f(B) = 1.5 and 2.5 do ; 
not show striking deviations from the observed values. 

Hovyever, it may he noted, for example^, for irradiation of 
line 6 with Vg = 2.2 Hz (Table V-2), that the calculated 
intensities for f(A)/f(B) = 1,5 are consistently higher than 
the observed ones for lines 1,2 and 3, and lower than the 
observed ones for lines 8,9 and 10, while the opposite is 
true for f(A)/f(B) = 2.5^ Such a feature is noticed for 
many different cases of irradiation. Tliis leads to the 
conclusion that f(A)/f(B) = 2.0 gives the best fit. The 
values for mechanism (ii), however, also agree q,uite well 
with the experiment. Similar agreement is obtained for 
irradiation of these and other lines in the spectrum for 
various strengths of irradiation, so long as the strength 
of irradiation is large enouigh to satisfy Bloch approximation, 
indicating that the relative importance of these two 
mechanisms for this molecule cannot be distinguished on the 
basis of this calculation, which depends only on the relative 
values of the elements of the relaxation matrix. 

The theoretical spectra for low’' values of Vg ^ 
however, depend on the absolute values of the relaxation 
terms. The observed and calculated relative intensities 
for irradiation of lines 6 and 7 with - 0.1 Hz are given 
in Table ('7-3), The calculated relative intensities for 
mechanism (i) are for f(A)/f(B) = 2.0, = 0.0 and three 

different values 0.1, 0.2 and 0,3 rad./sec. for 2f(B)T 
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These values are averages of ahout five speotra and are accurate approximately within 5 per cent 
of the values given. 
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(this was the oomiaoiL factor not included, in the calculation 

of relxation matrix elements for spectra mth large Vg) and 

—10 

those for mechaJiism (ii) are for t = 5.0 x 10 sec. and 
2.7 X 10“^^ sec."^ From the trend of the relative intensities 
given for mechanism (i) and hy mahing a similar comparison 
for irradiation of other lines, it v/as found that the spectra 
calculated with 2f(B)'i: =0,2 rad. /sec. show consistently 

good agreement. 

The results for mechanism (ii) with t = 5.0 x 10~^*^ 

- 0 

sec. also agree well with the experiment. The correlation 
time appropriate for this mechanism is that for rotational 
diffusion , for which a hard sphere model for the rotation of 
the molecule gives; 


- _ 4'%'na 

" “31cT 


(V-1) 


t These calculations vrere performed hy assuming tha,t only 
one off-diagonal element of % "between the states 
connected by irradiation, is non-zero (see Sec. IIG-2(a)). 
However for irradiation of lines in ' B ^ group in which each 
line has a close neighbour, more than one off-diagonal 
element of X- could become significant. To estimate the 
effect of this a calculation in the rotating :fram6., 
assuming two off-diagonal elements of % , between the 
ste^tes giving the close doublet (the doublet which is being 
irradiated), to be simultEaieniisly significant was per- 
formed (see Sec. IIC-2(b)). This calculation leads to 
values of X ’which are not significantly different from 

ixtx 

those with only one . off-diagonal element. 
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wiiere t] is the -visoosity and ''a'^ the hard sphere radius [?]« 
for the present s .ample p asJ 1.0 centipoise. At room tempe- 
rature, the value of 5.0 x 10“^^ sec. leads to a hard 

sphere radius of about 8.0 A° . This value of * a.' is consi- 
derably larger than the molecular dimensions, furthermore, 
the natural linev/idth of line 3j for example, given by 
R ^ //a; (full width at half height in Hz) becomes 0.65 Hz, 

for the above value of t , which is larger than any of the 
observed line?d.dths, including the inhomogeneity, for single 
resonance transitions. On the other hand, if, one-.; assumes a 
value of V’' of the order of the molecular dimensions 
(2 - 4 A^), the calculated intensities (see Table V-3i for 

T = 2.7 X 10”^^ sec. corresponding to ec^ual to 3.0 A ) 
do not agree with the experiment. These considerations show 
that mechanism (ii) is not a dominant relaxation process for 
this molecule. It should however be noted that the above 
argument is based on the estima.tes of correlation time for 
the hard sphere model of rotational diffusion, and therefore 
the values given should be considered within the limitations 
of this model. 

The effect of varying the correlation constant 0^ 
on the calculated spectrum for mechanism (i) is given in 
Table (V-4) along vmth the natiiral linewi dths ( R , /a ) of 

iX vX Ca# Uu 

the different transitions for different values of Also 

shown in the Table (Y-4) are the observed rela.tive intensities 
and the sum of relaxation and instability widths , obtained by 
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Table V-4 

Observed and oalonlated relative intensities (peal'i-beights) and 
linewidtiis of the double resonance treoisitions when line 6 is 
irradiated ’with an amplitude of 2,2 llz, The calculated inten- 
sities and linewidths are for mechanism (i) having C-p x, = 1.0, 
f(A)/f(B)=2.0 and 2f(B)T^=0.2 rad./sec. ^1^2 




Intensities • 



Linewidthst 




Galoulated 

Obse 

rve d^ 

Calculated 
relax, width 

S.l 

'll Obser 
. ved 


°iB 



V 



1.0 

0.5 

0.0 

Total 

Relax.+ 

instab. 

1.0 

0.5 

0.0 

1 

0.33 

0.29 

0.30 

0.30 

0.65 

0.29 

0.35 

0.31 

0.27 

2 

0,63 

0.61 

0.63 

0.64 

0.78 

0.42 

0.35 

0.30 

0.25 

'X 

0.42 

0.40 

0.41 

0.42 

0.83 

0.47 

0.34 

0.29 

0.24 

4 

0.54 

0.52 

0.52 

0.53 

0.54 

0.36 

0.24 

0.25 

0.25 

5 

1.00 

1.00 

1.00 

1.00 

0.60 

0.42 

0.13 

0.13 

0.13 

6 

- 

0.04 

0.04 

0.04 

- 

- 

0.24 

0.25 

0.25 

7 

0.46 

0.44 

0.44 

0.44 

0.62 

0.44 

0.26 

0.26 

0.25 

8 

0.20 

0.20 

0.20 

0.20 

0.50 

0.50 

0.15 

0,19 

0.24 

9 

0.28 

0.29 

0.29 

0.30 

0.42 

0.42 

0.17 

0,21 

0,25 

10 

0.09 

0.07 

0.07 

0.07 

— 

— 

0.19 

0.23 

0,27 

t 

Serial numbers, 

as in 

Table 

(Y-2). 






t Bull linewidths at half height in Hz, 

f^These values are averages over several spectra (5--8) and 
are accurate approximately within 5 per cent. 
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sulstr acting the inliomogeneity contribution from the obseirved 
widths, The intensities are not sufficiently sensitive to 
the value of C^, but the values of this 

constant quite considerably. It may be seen from the 
calculated v.alues of = 0.0, 0.5 and 1.0, 

that the effect of greater correlation between 1 .and B is to 
introduce an asymmetry in the natural linewidths in the 
spectrum, The observed relaxation plus- instability". Widths, 

however, do not show any systematic variations and since the 

the 

instability is assumed to contribute a same width to all the 
transitions a comparison of the observed and calculated widths 
shows tha.t is close to zero, Ihis estimate is only 
approximate as it is coupled with the estimates of inhomo- 
geneity contribution made from these linewidths. The effect 
of variation of 0^^ -□ is shovm in Table (V-5). This table 
shows tha.t neither the calculated double resonance relative 
intensities nor the calculated linev/idths are sufficiently 
sensitive to this parameter. However, as remarked earlier 


Og -Q has been taken as unity in the present case to maintain 
the complete equivalence of the two B protons. It may be 


noted that the values of correlation constants 



0 and 


-p =1.0 have been finally deduced in this work by 
"l^2 

experimental evidence which is not striking enough to make 


a sharp distinctiai between the different values. The values 


given may, therefore, be considered to be indica,tive of the 
degree of correlation, rather than to signify an accurate 
estimate of the same. 
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lable V-5 

Calculated relative intensities (peak-iieights) and 
relaxation widths of the double resonance transitions 
when line 6 is irradiated by an amplitude of 2,2 Hz. 
Ihese calculations are for mechanism (i) having 
f(A)/f(B) = 2.0, G ^ = 0,0, 2f(B)T = 0.2 rad./sec. 
and different valu^ of C-n -r, , 


Intensities Linev/idths 


S.H, 





^B^B2 


1.0 

0.5 

0.0 

1.0 

0.5 

0.0 

1 

0,30 

0.30 

0.31 

0,27 

0.24 

0.21 

2 

0.64 

0.65 

0.65 

0.25 

0.25 

0.24 

3 

0.42 

0.42 

0.42 

0.24 

0.22 

0.20 

4 

0.53 

0.53 

0.53 

0.25 

0.22 

0.19 

5 

1.00 

1.00 

1.00 

0.13 

0.16 

0.19 

6 

0.04 

0.04 

0.04 

0.25 

0.24 

0.22 

7 

0.44 

0.44 

0.44 

0.25 

0.22 

0.19 

8 

0.20 

0.20 

0.20 

0.24 

0.22 

0.20 

9 

0.30 

0.29 

0.29 

0.25 

0.25 

0.24 

10 

0.07 

0.07 

0.07 

0.27 

0.24 

0.21 
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The ■theoretical values of S”^? for the co.se of 
mechanism (i) ha'vlng f(A)/f(B) = 2.0, 2f(B)T =0.2 rad. /sec., 
= 0,0 and G^ ^ = 1.0 a.re plotted in juxtaposition vd.th 

the double resonance spectra, in Pigs. ('V’-3) - (V-7). There 
is good agreement between the theoretical and experimental 
spectra in all the cases. 

It shoul"d be pointed out that if the effect of 
magnetic field inhomogeneity is not properly tahen into 
account the values of the rela;;'Cation parameters determined 
and the possible distinction betv/een the theoretical calcu- 
lations for different mechanisms wo-uld be completely altered 
and may lead to erroneous results. For example, in the 
present case^if the magnetic field inhomogeneity is assumed 
to contribute the same vadth to the different double resonance 
transitions, the value f(A)/f(B) may have to be changed to 
nearly 4.0 and G^ to 0,5 to fit the observed spectra. 

It must be emphasized that it has been possible to 
conclude that mechanism (ii) is not an important relaxation 
process for the present case, only by combining the analyses 
of the results of double resonance spectra obtained for both 
strong and weah irradiation. The spectra for strong irra- 
diation contain, in general, more information in the sense 
that there are more transitions for ?irhich a fit can be 
obtained, but the final fit in this case depends only on the 
relat'ive values of the elements of the relaxation matrix. 

The spectra for weoh irradiation, on the other hand, contain 
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fewer parameters to obtain a. fit, but tiio calcnlation depends 
on the absolute values of the elements of the relaxation 
matrix. . ihis is an interesting distinotion betv/een these 
two categories of spectra and can be used with considerable 
advantage in the study of relaxation by double resonance 
methods. 

This analysis leads to the conclusion that the 
random isotropic fields at the site of nucleus A have a mean- 
square value twice tha,t at the sites of B nuclei. Ihis 
excludes the other spins in the molecule as important sources 
of these random fields, as nucletis A is farther off than the 
B nuclei from the remaining spins in the molecule. The 
proton relaxation thus appears to be caused primarily by 
intermole cular dipolar interactions. Since the sample is a 
20 per cent solution by weight in CCl^ (v/hich is nearly a 
saturated solution) each molecule of 2, 6-dibromoaniline is, 
on? an average, surrounded much more by CCl^ molecules than 
the molecules of the same kind (a 20 per cent solution by 
v;eight of 2, 6-dibromoaniline in CCl^ ha.s approximately seven 
GGl^ molecules per one 2, 6-dibromo aniline molecule). It 
was not feasible to look for a possible conceOi'fci^ation 
dependence of this relaxation since the range in which the 
concentration could be reduced, without seriously affecting 
the signal to noise ratio, would slightly increase the 
number of GGl^ molecules and would not significantly alter 
the physical situation in the immediate vicinity of a 
dibromo aniline molecule. .■ 
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VI - DOUBLE EESOIAJIGE SIUDY 01 PROIOI ESIAXAIM BT A IWQ - 
SPH SYSEEM (IE) COUPLED IQ A QUADRUPOLAR HUGLEUS 

VIA UTIRODUCIIOl 

In this chapter the resiilts of a proton-proton donhle 
resonance study of relazation in a two spin system (AB), the 
spectra of which exhibit broadening due to coupling with a 
neighbouring nucleus, are presented. The AB system is 
formed by the protons in 2-bromothi azole which is studied in 
the form of a neat liquid. The scalar coupling of the protons 
vd.th the nucleus modulated by the rapid quadrupolar rela- 
xation of the latter is an important source of proton relaxa- 
tion in this molecule. The analysis of double resonance 
expeilments performed by irradiating each of the single 
resonance transitions for different strengths of irradiation 
led to some detailed information on the scalar coupling 
interaction viz. the magnitudes and relative sign of proton- 
coupling constants, as well as the parameters describing 
the other mechanisms of relaxation in this molecule, in 
particular internal dipole-dipole interaction. Conventional 
relaxation time and linewidth measurements were also used to 
facilitate . the double resonance analysis. It was shown by this 
study that, in general, it is possible through relaxation 
studies by double resonance, to obtain the relative sign of 
the coupling constants of protons with q^uadrupolar nucleus; 

The spectra for '‘low’^ strengths of irradiation were analysed 
by the graphical method, mentioned in Sec. IIIC, 
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VIB BELATI7E SIGN OP SCALAR COUPLING CONSIiNIS 


When ■fciiere is more than one spin coupled to a quadru- 
polar nucleus, the relaxation matrix elements for the spin 
system contain terms arising from cross-products bstv/een 
various scsiar coupling interactions. These interactions are 
completely correlated iDecause the time-dependence is caused 
by the same quadrupolar nucleus and these cross-terms do not 
therefore vanish. For two spins both of which are coupled to 
the same q_uadrupolar nucleus the scalar coupling interaction 
may be written as (see Eq_. III-IO), 

I'A) + Jbj, I(B)>3,(t) 

The spectral densities (Eq. 11-30) will then contain terms 
of the type; 



1(A) + J 


BN 


I(B))*Sjj(t) 




(J^ 1(A) + JbnI(b)) 
(71-2) 


which depend^in general, oh the relative sign of and 

This suggests the possibility of determination of these rela- 
tive signs by studying the relaxation matrix elements. 


These relaxation matrix elements depend on the relative 
sign of and even in the case of single resonance. In 
single resonance spectra these matrix elements determine the 
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linewidths of the. different transitions. However , in cases 
in ?/hich ^ 

relaxation is contrih-uted only by the terns 27t[J^I^(A) + 
JBjjl 2 (B)]s|(t) of Eq. (VI-1) (see Sec. IIIG) and the widths 
of eanh of the transitions belonging in the zeroth order, . 
to 'a nucleus is (J^cos^0 + to the 'B'nucleus 

is 4 tc^ (J^ji-Gos^e + J^sin^0) where tan © = Jj^/( 

The single resonance spectra in such cases do not, therefore , 
depend on the relative sign of and Jgjj- 

Ihe double resonance spectra, however, depend on. a 
variety of relazation raatrix elements and are therefore 
expected to depend on the relative signs through the intensi- 
ties as well as linevd.dths of various transitions. This fact 
has been utilized in the present case to find the relative 
sign ■ J and Jb^* 

VIC EXPEEBIEHIAD PESUIIS 

1 . Single Resonance Spectrum 

Ihe single resonance proton spectrum of 2-bromo- 
thiazole is shown in Eig. (IV-5). Ihis is a typical four 
line AB spectrimi, having additional broadening due 'to sG.alar 
coupling of the protons with ni’crogen. The low field ('a" 
proton) doublet is assigned to the proton in position 4. 

This proton being nearer to nitrogen compared to the proton 
in position 5 (''B' proton) mil couple strongly with it, 
giving a broader doublet. This assignment is same as that 
in [44] . . - 
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The spectral ajialysis of an AB spin system is straight 
forward [2]. The analysis gives the values of the coupling 


constant = 3.55 + 0.05 Hz and the chemical shift 
= 18.8 + 0.1 Hz between the two protons. 


A B 


2. and Measurements 

Measurements of spin -lattice relaxation time and 

spin- spin relaxation time were made in order to derive 

some information on the magnitudes of scalar coupling of the 

protons with The magnetic field was made sufficiently 

inhomogeneous for the doublet of each proton to collapse into 

a broad line and and ^2 each such doublets were 

measured. ¥/as measured by the adiabatic rapid passage 

method [45]. Ihe result of one such experiment is shown in 

Big, (VI-1) and a plot be tween, log (M -M ) and time ''t"^ in 

^ z 

Big. (VI-2). The slope of this straight line is a measure of 
T^, The Values of T^ obtained by this method for the *A^ 
proton (low field broad doublet) and the S'"- proton (high 
field sharp doublet) are 9.1 + 0,5 sec. and 8,8 + 0,5 sec, 
respectively.. 

The spin- spin relaxation time Tg was measured for the 
'"B** proton by adiabatic rapid passage into the. center of a 
u-mode resonance with a large r,f, field (large compared to 
the inhomogeneity of the magnetic field), stopping the sweep 
at that point and monitoring the rate .of decay [45]. This 
decay is exponential with a time constant Tg and a typical 
curve obtained for the '^B^ doublet is shov/n in Big. (VI-3), 
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Fip; 


VI-1 . 


Eecovery . oux’ve for th- i 

wio adLiab.ntio rnn-,- -i 

method of T rr ' ■■‘■■'P-id pnC!iru;:o 

V '■ea.uro^.ont for the hj,,„ tun , 

doubipt of o T • ■’■•■'■‘■■-i-d 

psak-s are one aoeorfl , •»«-o-...rve 

■>^-oona apart-. 
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I'ig. VI-2. LogCM^ - M,^) versus time (t) graph for 
measurement of high field doublet of 
2-b romo thi az ol e . 




144 


Fig. YI-3. Saturation curve for the adialatic rapid 
pa,ssage method of Tg measurement for the 
high field doublet of 2~bromothiazole . 
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from which a value of I 2 obtained as 3«6 + 0.2 sec. The 
value of Tg ^" 0 ^ proton is comparatively very small 

due to its stronger coupling with and therefore could 

not he measured hy this experiment. 

nitrogen resonance was observed on a Varian 

V-4200 wideline n.m.r. spectrometer, operating at 2.695 MHz 

(field -^8.76 K.G.), using audio sweep modulation. Modulation 

frequency of 20 Hz and sweep amplitude of 1.0 gauss were used 

The spectrm was observed in the v-mode and the peak-to-peah 

linev/idth of the first derivative, after applying a modulation 

correction [46], was obtained as 0,98 gauss, which gives a 

fullwidth at half height of a lorentzian as 1,70 gauss. 

This leads to the relaXo,tion time t (T^ = Tp = t; ) for 

q^ r z q^ 

as 0,61 millisec. 

3* Do uble Resonance 

Some of the typical double resonance spectra obtained 
are shown in Figs. (VI-4)- (fI-6). The calibration for the 
strengths of irradiation was done in the manner described in 
Sec. VO. Figure (VI-4) shows the irradiation' of each one . 
of the single resonance lines by an amplitude of 0.135 Hz. 
'Low"^ irradiation amplitude spectra have been obtained for 
two more strengths viz. 0.076 and 0,043 Hz, for irradiation 
of each of. these lines. These spectra show significant 
intensity changes in the spectrum of the proton other than 
the one being irradiated. These changes are large for transi 
tions having energy levels in common with the pair being 
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Fig. YI--4. Frcquoncy swoop double resonaaioG spectra of 
2“bromothi azole, for '''low'' strengths of 
irradiation, obtained by irradiating (a) inner 
line of proton A, (line 3)5 (b) outer line of 
proton A, (lino 4), (c) inner line of proton 
B,(line 2) and ( d) outer line of proton B, 
(line 1). fhe strength of iriadiation in 
each c.ase is 0,135 Hz. 
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Pig. VI- 5. 


OlDserved and caloulalied dotillo retsomncc 
spectra for IrradLation of (a) outer lino of 
proton A (line 4) and (1) inner lino of proton 
A (line 3) in 2-1) romothi azole, with aniplifendo 
of irradiation 2.4 Hz. The calculated spectra 
are for internal dipole-dipole interaction 
and scalar coupling (mechanism (i)) with 
K = 4.0 and 0 = ^.3. 
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lig. 


YI“6, Observed and calculated double resonauce 

spectra of 2-bromothi azoic for irradiation 
of (a) outer line of proton A (line 4) vath 
strength of irradiation 1.6 Hz ond (b) outer 
line of proton B (line 1) with strength of 
irradiations 1.75 Hz. Iho cal-culated spectra 
are for internal dipole-dipolo interaction 
and scalar coupling (mechanism (i)) wi Hi 
K = 4-0 fpad == ^'5. 
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irradiated and can be q_ual*itatively imder stood on tbe basis 
of argument Sj, using an energy level diagram, similar to the 
ones described earlier (Sec. 7C-1). Ihe spectra v/ith * high^ 
strengths of irradiation show all the features of a double 
resonance spectrum (Figs. VI-5,6), The contribution to the 
linewidths to these double resonance transitions is also, in 
general, different for different transitions. Reliable 
estimates of this contribution for different transitions are 
obtained by using the plots of versus . These 

estimates are, hov;ever, approximate since the spin system , 
in the present case, is not wealcly coupled . 

Vlh MALTS IS 

1 • Eatiaates of Scalar Coupling Strengths 

The contribution of either dipole-dipole interaction 
or external random fields to T^ and T 2 , for very short corre- 
lation times, are equal (extreme narrowing), while this is 
usually not the case for the scalar coupling mechanism, since 
the corresponding correlation time is rather long in 
comparison [?]. Therefore^ wherever scalar coupling contri- 
butes to relaxation, a difference between the values of T^ 
and Tg gives information on the extent of this coupling. 
Estimates of coupling constants and of protons A and 
B respectively with may then be obtained by using Eqe. 
(111-49,50) and an assumption that the relaxation matrix 
element of the type represents (- 1 /T 2 )j^, where a — §' 
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is a ''A."' transition in the zeroth order^, This leads, for 
a coupled spin system AB under the assumption that only terms 
of the type I^S of in Eq, (VI-1) are significant, 

to: 


(1 - 1 ) 


2 1 




■ I ZL lg - ^t - U T^(J^ sin^e + cos^e) 

(VI-3) 


and 




= -- 4 - -^ ^ - ■'bh) J 

\ / \ iC 

where the subscripts A and B refer to the two protons and 
the subscript " sc ^ to scalar coupling. Equation (VI-4) 
gives the difference between the observed single resonance ' 
linewidths of A and B protons, if the contribution to these 
linewidths from other sources (namely, inhomogeneity of the 
magnetic field and other mechanisms of relaxation) is assumed 
to be the same. 

It may be noted that it is, in principle, possible to 
derive the values of the scalar coupling constants by double 
resonance data alone. However, estima.tes made on the basis 
of approximate determinations of and T2 sexve as a 
usefiil starting point for obtaining fit with the double 
resonance experiments, and v/ould facilitate the determina- 
tion of the parameters of the other mechanisms. 
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Using the experimentally determined values of 
and (l2^B ^®spectively 8,8 sec. and 5*6 sec. and = 0.61 
nillisec, in Eq_. (VI-3)» and equating the observed linewidth 
difference between the A and B proton resonances (= 0.35 Hz) 
to Eq. (VI-4), gives: = 8.8 Hz and Igp = 3.1 Hz. This 

value of compares favourably with a similar estimate of 
10.8 Hz by Kintzinger and Lehn for thiazole [37]. 


For these VEjlues of and and the observed % 

the condition mentioned in Sec. IIIC, for complete washing 
out of the scalar coupling 1) is satisfied. Therefore_,for 

the present sample the treatment of scalar coupling only as 
a mechanism for relaxation is valid. 


2. Other Details 

The proton-proton distance in this molecule enters 
the internal dipole-dipole interantion calculation. Informa- 
tion is available on bond angles and bond lengths in thisizole, 
through studies of 0 ^-H coupling constants [47]. Assuming 
that the bromine substitution does not significantly alter 
the geometry of the molecule, the distance between the tv/o 
protons in 2-bromothi azole is obtaine-d ars 1.66 A*^ . 

Theoretical single and double resonance spectra have 
been calculated by considering (i) internal dipole-dipole 
interaction and scalar coupling and (ii) external isotropic 
random fields and scalar coupling as possible alternative 
mechanisms for relaxation. ■ For double resonance spectra with 
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''high'' strengths of irradiation, for which jYH 2 | | {~aa^ 

Bloch approximation has been used (see Sec. IIG). Since 
under Bloch approximation the calculated spectra depend only 
on the relative value of the relaxation matrix elements, 
these spectra have been oalciiLated by varying a parameter K, 
which determines the relative contribution of dipole- dipole 
<and scalar coupling in mechanism (i) and random field and 
scalar coupling in mechanism (ii). Bor random fields, however, 
there are two more parameter namely the ratio of mean 
square fields f(A)/f(B) at the sites of the two protons and 
0 the correlation constant between these fields. The ratio 
f(A)/f(B) was varied in the calculations to obtain a fit with 
experiment, but C was fixed from the information available 
from the observed single resonance linewidths. Ihe observed 
single resonance linevddths are more for the two inner lines 
in the AB quartet* in this ca.se, than the outer lines, indi- 
cating tha,t G is close to unity [ref. 7, p. 509], and was 
tahen to be 1. It may f'urther be remarked in this connection 
that while the internal dipole-dipole also contributes greater 
Viddths to the inner lines compared to outer ones, in the AB 
quartet, the scalar coupling contributes equal widths to the 
two lines of a proton resonance in it. 

Bor the double , resonance spectra with low strengths 
of irradiation, the observed and calculated straight line 
plots (Eq, 11-80) were drawn and compared. 
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YIE Discussion 

Table (VI-1) sliows the observed ajid calculated double 
resonance peak heights, the observed single and double reso- 
nance linewidths and calculated single and double resonance 
relaxation widths for irradiation of outer line of ' reso- 
nance by anplitudes of irradiation V 2 = 1.6 and 2,4 Hz and 

\ 4 

for irradiation of outer line of B resonance by an ampli- 
tude Vg = 1.75 Hz, A satisfactory fit for observed peak: 
heights and linewidths is obtained for mechanism (i) for K = 
4,0 and = 0,5 and for mechanism (ii) for K = 4.0, 

= 0,3 and f(A)/f(B) = 0,5. The parameter K is equal 
to [ mechanism (i) and equal 

to [ (16m^ J^Tg/3)/(2f(B)i:^)]® for mechanism (ii), A variation 
in does not seriously affect the calculated peah- 

heights, but alters the calculated linewidths. The 

versus A- plots for approximate estimation of inhomogeneity 
J 

contribution to vailous double resonance transitions were 
drawn and one such plot for irradiation of 'a" proton by v^ = 
2.4 Hz is shov/n in Pig, (YI-7). This plot indie ates, for 
irradiation of the outer line of 'A'" resonance which corres- 
ponds to A^ =-1,94 Hz, approximately 40 per cent increase in 
the recoverable part of single resonance inhomogeneity 
contribution (the instability contribution is not recoverable) 

f Though this quantity has been estimated from T^ and T^ 

measurements, the accuracy of these measurements is limited 
in as much as these measurements were made for overlapping 
doublets. 
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Mg. 


Plots of ca3-ctilated double rosouauee rrociii<.:u-- 
oies of 'B’' transitions versus ''freciiunio;/- 
offsef^of irradiation of ''a'^ nuclei, ihesc 
plots are oolculated foi t.no KS> spin syoteri 
for Vg = 2.4 Hz. Irradiation of line 1,2,3 
and 4 oorre spend respec’bively to = 20.74 > 

17.19 , 1*61 “1.94 Kz. 
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for lines 1 and 2 of the double resonance speotrutn and a 
corresponding decrease in lines 3 and 4. If the recoverable 
part of the inhomogeneity is assumed to contribute a full- 
mdth at half height 0.15 Hz to single resonance transi- 
tions (see Sec. 70-2), the observed and calculated single 
and double resonance mdths fit best for = 0.3* A 

detailed comparison of the double resonance spectra obtained 
for different frequencies and amplitudes of irradiation 

-juiajj led to a satisfactory fit for the above mentioned 
parameters for mechanisms (i) and (ii)- 

The effect of a change of sign of -with respect 
to is also shorn in Table (VI-1), The single resonance 
line-widths are not changed as sho-vm in Sec, (VIB). However, 
the double resonance linewidths shov/ a marked difference 
between the two alternative cases of same sign and opposite 
sign. The linewidths for the case of opposite sign are widely 
different from the observed ones, leading to a conclusion that 
both and are of the same sign. 

The spectra for * low ^amplitudes of irradiation have 

been calculated for the above two mechanisms with parameters 

K and same as]fw1iioh led to a satisfactory fit for 

' high' amplitude of irradiation case, and by varying the 

absolute values of the relaxation matrix elements. This 

amounts to, for a given value of t , a simultaneous variation 

Q, 

in and for mechanism (i) and in 2f(B)'T:^, and 

for mechanism (ii) such that K and are unchanged. 

A fit for '‘low^ and 'high' amplitude of irradiation spectra; 
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then determines the values of these parameters as obtained 

by double resonance analysis. However, some information on 

available from relaxation time measurements 

and therefore the absolute values of the relaxation matrix 

elements need be varied only over a small range. Table 

(YI-2) gives the observed and calculated values of R, 

Eq_, ( 11-79 )j for various transitions in the double resonance 

spectra for ' low ^ irradiation of A group of lines. The 

calculated values axe for K = 4.0; = 0.3 and 

different absolute values of relaxation matrix elements. A 

satisfactory fit with the observed values of R is obtained 

for mechanism (i) with t = 7.0 x 10“ sec. For mechanism 

(ii) no satisfactory fit is obtained. Also shorn. in Fig, 

(VI-8) are the plots of observed and calculated ratios 

(R+1)/(R-1) between various transitions versus l/(v|), for 

irradiation of the ^A^ group of lines with ''lo^’v^ amplitudes 

of irradiation (see Eq_. 11-80). While the observed ratios 

are indicated by va,rious points on the graph, the calculated 

ratios are shown by straight lines. The fit between the 

observed points and the calculated ra,tios for mechanism (i) 

(continuous lines), for K = 4.0, = 0.3 and = 

—12 

7.0 X 10 sec, is satisfactory. These parameters give 
'^AM ^ *^31 *^2.9 Hz for = 0.61 millisec. For 

mechanism (ii), however, no satisfactory Fit could be 
obtained for any plausible values of the parameters, and the 
plots in Fig. (VI-8) for this mechanism (broken lines), 
for K = 4.0; = 0,3 and 2f(B)r^ = 0.19 rad./sec. 
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Talole VI-2 

^^observed and calculated values of R (sec Bq. 11-79) for. 
low" irradiation of 'A' group of lines, 'Ihe calculated values 
for both mechanisms fi) and (ii) are for K = 4,0 and = 

0,3. Bor mechanism (ii) the additional parameters are f(AT7f(B) 
= 0,5 and G =, 1.0, 


Irrad- 

iation 

at 


Vo Ratio 

2 - Oboe 

in Hz T rvod 

1 Y) n p. 


mechanism (i) 
_12 

Vq X 10 sac 

11.0 7.0 4,7 


mechanism (ii) 
2f(B)v^ 

0.29 0.19 0.13 


2/1 0.76 0.88 0,77 0.64 0.94 0.39 0.81 

0.043 

3/1 0,89 0.94 0.88 0.82 0.93 0.87 0.77 

2/1 0.61 0 i 72 0.59 0.48 0.86 0.76 0.68 

line 4 0,076 

3/1 0.79 0.86 0.79 0.73 0.83 0.73 0.63 

2/1 0.50 0.54 0.45 0.40 0.73 0.64 0.59 

0.135 

3/1 0.71 0.76 0.71 0.69 0.69 0.59 0.53 

2/1 1.33 ia8 1.36 1.67 1.06 1.12 1.23 

0.043 

4/1 1.24 1.10 1.20 1.38 0.98 0,96 0.92 


2/1 1.89 1.47 1.84 2.30 1.16 1.28 1.43 

line 3 0,076 

4/1 1.48 1,27 1.48 1.74 0.94 0.90 0.84 


2/1 2.54 2.04 2.50 2,88 1.33 1.48 1.60 

4/1 1.89-1.59 1,86 2.07 0.88 0.81 0,78 


0.135 
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Fig. VI-8, Plots of observed ond caloulo-ted ratios (R+l)/ 
(R-1) versus 1/^1 for^ low^ strengths of irradia- 
tion, for 2-bromothiaKOlG . 'i'he point; s indioato 
the observed ratios onci the straight" linen the 
theoretical ratios. The solid linos are for 

mechanism (i) with K = 4.0, “ 0*3 

_1 o 

“ 7-0 X 10 sec. the broken linen are 

for mechrnism (ii) with K = 4.0, ^ 0.';:), 

2f(B)i;^ = 0.19 rad. /sec., f(,A)/f(B) = 0.5 and 
^iU3 ” 1»0. For irradiation of lino 3 A A 
represents, in the defirdtion of 11 , tluj ratio 
of lino 4 to lino 1 and SB the ratio of lint; 

2 to line 1, For irradiation of lino 4 OC 
reprrasonts the ratio of lino 2 to line 1 and 
Ah the ratio of line 3 to line; 1. For each 
place of irradiation three nti’engths 0.043, 
0,076 and 0.135 Ha were used. 



{R + l)/(R- 1) 
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i ^or ^ 2.9 Hz) show large deiriations 

from experimental values. Ihese considerations show that 
the external random fields along xvitii scalar coupling do not 
account for the ohserved changes in the spectra with low'^ 
amplitudes of irradiation and hence mechanism (ii) is ruled 
out. The internal dipole-dipole interaction along vri.th 
the scalar coupling accounts for all the features of the 
single and double resonance spectra, and is ,therefore, the 
dominant mechanism for relaxation in tais molecule. As the 
relaxation contribution of internal dipole-dipole interaction 
is strongly dependent on the distance between the interac- 
ting spins (r"^), this result is not su3:?prising in view of 
the small proton-prpton distance (1,66 A®), in this molecule, 
(compg,red to 2.48 betv/een the ortho protons in benzene). 
Ihe viscosity of 2-bromothi azole was measured at room tempe- 
rature and a value ^ 2,0 centipoise obtained. This leads 
to a hard sphere radius (see Eq. V-1) of the order of 2 A^ , 
which fits well v,ath the molecular dimensions in this case. 

The theoretical spectra for mechanism (i) having 
K = 4.0 and - 0*3 are plotted in Juxtaposition with 

the double resonance spectra in Figs. (VI-5 and 6). These 
show good agreement between the observed and calculated 
spectra. 

Table (VI-3) gives the contribution of internal 
dipole-dipole interaction and scalar coupling to the single 
resonance relaxation matrix elements of the type 


and 
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Table VI-3 

Contributions of internal dip ole- dip ole interaction 
and scalar coupling to various single resonance 
relaxation matrix elements (in rad, /sec*) for meciia- 
nism (i) having K = 4.0, = 0.3 and = 

7.0 X 10"^^ sec. 


Matrix 

elcnentt 

Bipole- 

dipole 

Scalar 

coupling 

Total 

12 12 

-0.174 

-0.149 

-0.323 

13 13 

-0.149 

-1.508 

-1.657 

14 14 

-0.171 

-2.565 

-2.737 

2 3 2 3 

-0.076 

-0.731 

-0.808 

2 4 2 4 

-0.174 

-1.508 

-1.681 

3 4 3 4 

-0.149 

-0.149 

-0*298 

1111 

0.171 

0.000 

0.171 

112 2 

0.035 

0.000 

0.035 

113 3 

0.021 

0.000 

0.021 

1 1 4.4 

0.114 

0.000 

0.114 

2 2 2 2 

-0.086 

-0.013 

, -0.099 

2 2 3 3 

0.018 

0.013 

0.031 

2 2 .4 4 

0.034 

0.000 

0.034 

3 3 3 3 

-0.065 

-0.013 

-0.078 

3 3 4 4 

0.023 

0.000 

0,023 

4 4 4 4 

-0.171 

0.000 

-0.171 


The numbering of the states is such that for the 
iB spin system, in the zeroth order these states 
are given, in the increasing orde^; by, aa, ap, pa, 
and pp. 
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for the parameters for vihlch a satisfactory fit of 
the double resonamoe spectra has been obtained, This table 
shows the relative contribution of these two mechanisins to 
the relaxation in this spin system. 

It is concluded from this study that the internal 
dipole-dipole interaction and the scalar coupling of the 
protons with nitrogen, v/hich is relaxing rapidly due to 
quadrupole relaxation, in 2-bronotlii azole, accountsfor proton 
relaxation in this molecule and that the external random 
fields do not contribute significantly to the relaxation 
processes. Various parameters describing the internal 
dipole-dipole interaction (r 7..0 x 10 sec.), and the 
scalar coupling mechanism ^ 9.7 Hz, ^ 2.9 Hz and 

of the same sign) have also been obtained. 
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SDM&IAEY MD GOxIOlUSIONS 


Relaxation processes in oonpled spin systems have 
been studied by steady-state hoiaonuclear double resonance 
experiments. The systematics of the double resonance method 
and details of the relaxation processes have been examined. 
Double resonance experiments have been performed on a coupled 
s3n3metriGal three spin system and on an iB spin system. 

It is found that it is useful to study double 
resonance spectra obtained over a wide range of strengths of 
irradiation, as the solution of the density-matrix equations 
obtained for conditions of high and Iot; irradiation 
amplitudes possess some distinctive features v/hich may then 
be used with considerable advantage to distinguish between 
different relaxation mechanisms. A detailed consideration 
of the inhomogeneity of the magnetic field is shown to be 
of crucial importance for the analysis of the relaxation 
effects and, a method of including it in the density matrix 
calculations of the double resonance spectra is' suggested, 

Ihe inhomogeneity of the magnetic field is found to contri- 
bute a ftill -width at half-height of the order of 0,15 Hz to 
single resonance transitions^ in the present experiments, which 
was almost wholly recoverable in double resonance spectra 
under suitable experimental conditions. 

Isotropic random field model, convenient for descrip- 
tion of relaxation effects in such studies^ is shown to be 
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a plausible model for various iut enrol ecular and intramole- 
cular interactions under conditions normally realised in 
experiments. Whenever such a model is applicable the cons- 
tants describing the correlation between the random fields 
at the sites of different spins have been discussed. A 
.definition of this constant ,v/hich includes the possibility 
of completely correlated fields with different mean square 
values, is used. 

The relaxation of the AB^ spin system of the ring 
protons of 2,6-dibromoaniline, -which was studied as a 20 per 
cent solution in CCl^, is found to be predominantly due to 
sxtemal random fields having mean square value tmce at 
the site of ' pro'ijon compared to either of the two '' B ^ 
protons. The random fields at the tv;-o 'b'' protons are found 
to be highly correlated, v/hile those at 'A'' and any of the 
""B^ protons being almost unoorrelated, Intermolecular dipolar 
interaction between the protons and the chlorine nuclei of 
the solvent appears to be the source for these random fields. 
Ihe effect of intramolecular dipolar interactions between the . 
various protons, ¥d-thin the spin system, is small compared to 
that due to the intermolecular random fields, 

Ihe relaxation of the AB spin system of the protons 

of 2-bromothiazole, studied as a neat liquid, is, on the other 

hand, found predominantly by the intramolecular dipolar 

interactio^n between the two protons and by their scalar 

14 - ^ 

coupling with the nearby h which intum is relaxing fast 
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due... to quadrupolar relazation. liie intrajnolecular dipolar 


interaction is significant in this molecule, because the 

I ihai 

proton-proton distance is small comparedtOAin molecules like 
benzene. External random fields do not seem to contribute 
significantly in this molecule. Double resonance experiments 
performed over a mde range of strengths of irradiation and 
at various frequencies of irradiation led to the determina- 
tion of the magnitudCLS.. and relative sign of proton-U^^ 
scalar coupling constants. Ihe various parameters describing 


.-1 9 

the proton relazation are, x = 7.0 z 10 sec., J c 

^ G * Ml ' 

Hz and ::r'. 2.9 Hz (for x^ = 0.61 millisec)? and J 


9.7 


BH 

have the same sign. 


9 


BH 


It may be noted that the presence of bromine does 
not modulate proton relazation in both the above cases. Ihe 
random fields arising due to the scal.ar coupling of protons 
and bromine are negligible presumably due to a fast quadrupole 
relazation of bromine nuclei which possess a large quadrupole 
moment. 

Further studies of this hind can be aimed at (i) 
studying the relative importance of dipole-dipole interaction 
and eztemal random fields as a function of distance between 
the interacting spins, (ii) studies by double resonance in 
which scalar coupling. is only partly '' washed-out'^, (iii) relaxa- 
tion effects in double resonance as a function of external 
parameters like temperature, external magnetic field etc., 
for determining the relative inportance of various 
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interactions, in particular, that of the intramolecular and 
intemolecular dipole-dipole interactions, spin-rotation 
interaction, anisotropic chemical shift (if any). liirther 
insight into the relaxation processes can also be gained by 
attemptirjg to interpret the different values of the mean 
square random fields and their mutual correlations in terms 
of specific models for the dynamical behaviour of the 
liquid. 
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APPENDIX I 


LISTING OF COIIPUTER PROGRAI'/I 



BTTRODUCTIOH 


The listing of a typical computer program used in the 
density-matrix analysis of the double resonance spectra, is 
given in the following. The program is written in Fortran IT 
and run on IBM 7044. This program in its present form is 
applicable upto three spin systems, but can easily be extended 
to larger spin systems. Further, this program calculates dotble 
resonance spectra for high values of V 2 , for which Bloch 
approximation holds. For *low* strengths of irradiation, the 
laboratory coordinate calcxfLations have been performed by 
incorporating minor modifications to include the off-diagonal 
matrix element of pC . 

In the first part of this program the spin hamiltonian 
matrix is formed in the spin product basis, A matrix dia- 

gonalization routine is used to obtain the eigenvalues."^ This 
routine also gives the transformation coefficients which are 
then used to transform various quantities from spin product 
basis to a,a'^.,. basis, in which is diagonal. 

In the second part of the program, various relaxation 
matrix elements (R) are calculated by using a subroutine named 
RFDCACi. The relaxation matrix elements can be calculated for 
any degree of 'correlation. The simultaneous equations for 
matrix elements of % are then formed and solved by using a 
determinant evaluation subroutine named liBTER. In the last 
part the relative intensities of the double resonaace transi- 
tions including relaxation effects are calculated. 

+ This routine and a part of the remaining program are similar 
to those used by the magnetic resonance group at the Depart- 
ment of Chemistry, Harvard University, Cambridge, Mass., 
U.S.A. during 1961-65. 
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SJOB PHG'JC'7, TIHEOOS, PAGESC30, MAME ANIL 
STB JOB 

SIdFTC MAI “4 •■40DECK 

C TF4S pnoGRAM CAlCjlATEo DOUSlE RESONANCE SPECTRA USING EXTREME 
C NARROWING ASSUMPTION AND BLOCH APPROX I ‘^‘AT I ON . 

C r: 0 SPlN IS THE order qF THE SPIN SYSTEM CONSIDERED. NASSYM IS THE 
C NUMBER OF ASSYMMETRIC STATES. CHART IS THE SPEED OF CHART PAPER 
C PER C.PoS. And gives direct CALIoERATION of THE SPECTRA. 

C THE F ARE COJPONANTS OF LATTICE FUNCTIONS IN THE IRREDUCIBIlE FORM 
C THESE may be EXPRESSED WITHIN A PROPORTIONALITY CONSTANT 
C ;'C.;TAX IS THE NUMBER OF COMPONANTS OF LATTICE OPERATORS. 

C MS STORES THE MATRIX ELcmEnTS OF THE RclAxATIhn HAMILTON IAN IN. THE 

c product basis, fed as Input data, the chemical shift (A-bi, 

C Ti E coupling constant R, THE AMPLITUDE OF IRRADIATION VA=VB AMD THE 
c Frequency oi irradiation ia=frEq.of a resonance- frEq. of irrao. 

C AND o^FrEq. oF B rESOnAnCE - FrEq, oF IrrAD.) ArE INPUT PARAMETERS 
C Cl, C2» and C3 are THE DEGREE OF CORRELATION AND ArE TO BE FED AS 
C H PUT PARAMETERS. 

C F' M 0 R FOR BLOCH APPROXIMATION 
C Ni'iDR three SPIN 1/2 

DxMENSION ABC{3,8,15i ,F(15) ,As(8,8,15) ,ROOT{10) ,AA(10»10) 

dimension X{8,8),Y(a,8),E(8,8),H(8) 

dimension UdO.lO), V(8) ,SIGMA(10,10) ,REPVECdO) 

dimension RAAB3( 10,10) ,RABAB(S,8) ,RAA3C(8,8,8) 

common NOSP I N , NQMAX , F 

C 0 M CM / A N I L / A B C , C 1 » C 2 , C 3 

COMr-.ON / K UMAR / NS YM 

STW0=SQRT{2.C ) 

READ1O07, CHART 
PRiNTlOlO, CHART 

A RFAD1009,NOSPIN, NASSYM 
NDlM = 2*^:-N0LniN 
ISYM=NDIM-NASSYiM 
PRINT1006 

7 READ1005 ■ 

PRIMT1005 
READ1009, NQMAX 
D012')NQ=1, NQMAX 
DO 12C I=1,NDIH 
DO 120 J=1,NDIH 

120 AS( I ,J,NO)=0.0 
DOl2iNQ=l ,NQMAX 

11*^ READi030, I ,J»ASS 

IF{ I ,EQ.O)GO TO 121 
ASI I ♦J,NQ)=ASS 
PRINT1031,I,J,NQ,ASS 
GO T0119 

121 CONTINUE 

8 READ1Q07, IFd ) »I = 1»NGMAX) 

PRINT1007,(Fd ) »I = l*NQiMAX) 

c xd.j) Are THE matrix Elements of the spin hamiltonian in the spin 

C TROOUCT BAST’S ^ 

C For DIFFErEi.T spin systems the following gards have to be changed 

C THE present SET OF CARDS ARE FOR AN AB2' TYPE OF SPIN SYSTEM. 

C FHE ANTISYMMETRIC STATES, IF ANY, SHOULD BE THE lAST PAIR OF STATES. 



1 RF'D 1007,Ai,B,R»VA»VB 
PR1NT1006 

’RINT10C0,a»B»R»VA,VB 
C _ COMPUTE X MATRIX 
^ DO 

D05.t = l,NDIM 
5 xn»J)=o. 

X( T,l) = ( a+2»-»B+R)/2. 

X ( 2 > 2 ) = { -a+ 2 o^«-B-R ) /2 , 

X{3,3)=A/2, 

X{A,A)=-a/2. 

X { 5 ,5 ) = ( A-2o*8-R 5 /2> 

X(6»6)=(-A-2b*B+R)/2. 

X{7,7)=A/2o 
X( 8 »8 )=-A/2» 

X( l,3)=VB/STWO 
X(l,2)=VA/2, 

X (2,4)=\/b/STWO 
X{3»5)=VB/STWO 
X(A,6I=VB/STWO 
X(3,A)=VA/2. 

X(4,5)=R/STWO 
X ( 2 » 3 3 =R / STWO 
X(5,6)=Va/2. 

X(7,fa}=VA/2. 

DO 3 I=1»NDIM 
DO 3 J=I,NOIM 
3 X(J,l)=X{ I,J) 

C 

c i-iATRIX DIAGCNALIZATION ROUTIiME 

C arrays USED X(I»J)» E { I ,J ) ( EIGNVECTORS) , D I MENS I ON-NDI M 

C FIND LARGEST OFF DIAGONAL X 

C 

25u DO 233 I=1,NDIM 
DO 253 J=1,NDIM 
IF(I-J) 251,252,251 

251 E(I,J)=0, 

GO TO 253 

252 E(I,J}=1,0 

253 CONTINUE 

2 00 LA=1 
LB = 2 

I3X = ABS{X(2,1) 3 
DO 202 J=3,NDIM 
JA=U-1 

DO L02 1 = 1, JA 
D=3X“ABS(X( I»J) 3 
1F(D3 201,202,202 

201 0X=ABS{Xf I ,J) 3 

la=i 

lb=j 

20? CONTINUE 

IF(BX-0. 000013 ' 320*320,300 
C COMPUTE ROTATION 
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30'* BX = X(LA»LE)*X«LA»L6) 

D=X(LA>LA)-X{LB»LB) 

R = SQRT(D*D+4,-;i-3X) 

A=SGRT ( APS { ( R + D ) / ( 2*C--^-R ) ) ) 

IF(C.707-A) 302,302,301 
3nl 3=n.C-A 

A=SQRTa»0-B*8) 

GO TO 3n3 

3r'2 B = -SQRT ( 1.0-A-»A ) 

3k.3 IF(D/X(L/ ,LB) ) 304,310,310 
304 b=0.0-B 

c orthogonal rotation of matrix 

310 DO 313 J=1,NDIH 

IF(LA-J) 311,313,311 
313 IF{L8-J) 312,313,312 

312 Xn.A,j)=A*X(J,LA)-B*X(J,LB3 
X ( lb » J ) =B*x ( j,LA )+a^;-x ( j, lb ) 

313 CONTINUE 

DO 314 J=1,NDIM 
X( J,LA)=X{LA,J) 

31-+ X( J,LB)=X(LB,J) 

D=X(LA,LA)+X(LB,LB) 

X{ LA,LA) =A*A*X (LA,LA)+B*B*X{LB,LB)-2,*A*B*X{LA,LB) 
X(LE,LB)=D-X(1 A,LA) 

X{LA,L6)=0.0 

X(LB,LA3=0»0 

C Ep-^ER ROTATION IN E MATRIX 

DO 315 J=1,NDIM 
3X=E( J,LA) 

E ( J,LA)=A*BX-8*E £ J,LB ) 

313 E£ J,LB) = B^f-BX+A-»-E£ J,L8) 

GO TO 200 

C TEST DIAGONALIZED MATRIX 

320 EX -0,00004 
JA = NDI 1-1 

DO 323 1 = 1, JA 
JB=I+1 

DO 323 J=JB,NDIM 
IF{X£ I,J) ) 321,323,321 

321 R=ABS(X(I ,J)/(X(I,I)-X(J,J> ) ) 

IF(R-BX} 323,323,322 | 

32? BX=R 
LA=I 
L8 = J 

323 CONTINUE 

IF(BX-0. 00005} 450,450,300 
C CALC TRANSITION INTENSITIES AND FREQUENCIES 

4^0 DO 451 I=1,NDIM 
H(I)=X{I,I) 

DO 451 J=1,N0IM 
451 XCI ,J)=0.0 
C 

c THE Following are matrix elements of k + i in spin product basis. 
C these CARDS have TO BE CHANGED FOR DIFFERENT SPIN SYSTEtMS. 
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c 

X ( 1 , 2 ) = 1 e 0 

X(1,3)=STW0 

X{2,4)=STW0 

X(ji54)=lo0 

X(3,5)=STWO 
X { ^ » 6 ) = S T .j'O 
X ( 5 » 6 ! = 1 
X(7,8)=l, 

452 DO 453 I = l,iMDlM 
DO 453 J=l,NDi:‘1 
Y( I , J)=0,0 

DO 453 LA=1>NDIM 
DO 453 l8=3»NDIM 

453 Y( i ,J)=Y( I ,.])+ E(LA,I )*X{LA,LB)i^E(LB,J) 

:’RINT 1001 

DO 456 I=1,NDIM 
PRINT 1002,I,H(I) 

DO 456 J=1,NDIH 
R=Hl I )-H( J) 

YS0='{ I *J)*Y( I »J) 

plot=r>chart 

456 PRINT 1003, J,R ,YSQ,Y( I , J ) ,E{I,J3,PL0T 
C CALCULATION OF REDFIELD ELEiYEFiTS 

C TRANSFORM RELAXATION MATRIX 

D0130nq= 1 ,NQHAX 
DO 13o I=1,NDIM 
DO 13 o J=i,NDlM 
ABC( I ,J,NQ)=0.0 
D0130(< = 1 ,NDIM 
DO 13 o L=1,NDIM 

13U ABC{ I ,J,NQ)=ABC(I,J,HO)+E(i<»I )*AS(K»L»N03^-E{L,J) 

C FORM SIGMA MATRIX AND TRANSFORM IT 
00515 1=1, NDIM 
D0515J=1 ,NDIM 
5x5 X( I » J)=0.0 
C 

C THE FOLLOWING GIVE -F < Z ) =-SU!i ( 1 ( Z ) ) IN THE SPIN PRODUCT BASIS. 

C THESE Are THUS MATRIX ELEMENTS OF SIGMA ZERO IN SPIN PRODUCT BASIS. 
C THESE Cards HAVE To BE CHANGED FOR DIFFERENT SPIN SYSTEMS. 

C 

X(l,l)=-1.5 
X<2,2)=-0.5 
X{3,3)=-0.5 
X(4»4)= 0.5 
. X(5,5}= 0.5 
X( 6,6)= 1,5 
X(7,7)=-0,5 
X(8,8)= 0.5 

C TRANSFORM SIGMA I.ATRIX 
PRINT1020 
D05i7l=i ,NDIM 



I , J) = 0.0 
,>05i6L = l ,.X‘DrM 
DOS j.&i<.= l ,NDIM 

5 io oIG‘»iA( I , J ) = SIGi-'!A ( I » J ) +£ ( K. , 1) ( K » L) *E f L » J ) 

517 PRINT1003, (SI6i'.1A ( I , J ) , J= 1 , MD IM ) 
pRMTioas 

c rpad CORRELATON 
6 READ 1.>02,NC0RLN 
D0900 IJK=1,NC0RLM 
C Cl, C2 and C3 are CORRELATIONS. 

9 READloo7,Cl ,C2 ,C3 
A5P PRINT 1006 

PRINT1024,C1,C2,C3 
PR I NT 1006 

C CALCULATION OF R ( n J J ) 

DOSlC 1=1 ,NSYM 
D0509j=i , nsYM 
CALL RFDGAL( I ,I ,J,J,RFD) 

HAA6B( J, I )=RFD 
509 RAABBC I ,J)=RFD 
5:' 9 CONTINUE 

D0512j = l ,NSYi'1 
D0512K=J ,i4SYN 

IF(JoEQ.K) GO TO 512 

00511 1=1 ,NSYM 

CALL RFDCALU , I , J,K,RfD) 

RAABC( I , J,K )=RFD 

511 R/ ABC ( I »K»J)=RFD 

512 CONTINUE 
D0514i=l ,NDIN 
D0513j=l ,NDIM 

IF( I.EQoJ) GO TO 513 
CALL RFDCAL( I , J, I ,J,RFD) 

RABAB( I , J ) =RFD 
RABAB( J,I )=RFD 

513 CONTINUE 
11=1+1 

51 A. PRINT1015 , ( (RA8AB( I » J ) , I , J , I , J ) , J= 1 1 ,nDI M ) 
C SOLU^OION OF simultaneous EQUATIONS 
D0S20I = 1 ,NsYM 
n n=o. 

D0520j=l ,NSYM 
D052GK=J,NSYM 
IF{J,EQ.K)60 TO 520 

n=V( I )+2.0*SIGMA( J,K}-»RAABC( I ,U,K) 

520 CONTINUE 
PRiNTlllS 
D052Al = l ,NSYM 

52A PRINT1019, {RAABB( I »J) »J = 1,N5YM) ,V(I ) 
PRINT1006 

CALL DETER(RAABB,DET) 
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PRiNTloI i »DET 

C CHECK, FOR LINEAR DEPENDEWCE 
00530K=1 ,NSYM 
D052d I = 1 , KiSYM 
D0525J=1 jNsyM 
rfpveC(I)=V(I) 

52 5 U( I 9 J )=RaaBB{ I » J 5 
D0 5 26j=l jNEYiM 
526 U(K»J)=l.o' 

REPVFCIK ) =0,0 
CALL DETTR (U,DET } 

PRlNTlo2i,K,DET 
CHFCK=A3S( DFT ) 

IF ( CHECK-0.0001 > 530 9 540,54'' 

530 CONTINUE 

54r IF ( CHECK, LT,C, 0001 ) GO TO 900 
DC 5b'.-K = l 9NSYH 
D0571 J=1 9MSYM 

D0571 1=1 jNSYM 
AA ( I 9 J) =U{ ! ,J ) 

IF(J-K) 571,570,571 
57u AA( I »K)=R£PVEC( I ) 

371 CONTINUE 

CALL deter (AA, DEM) 

R00T{K)=DEM/DrT 
PR INTI 006 

580 PRlNTl0l2,K,DEM,K,R00T(K) 

C CHEr'K the ROOTS 
PR I NT 1006 
^RINT1016 
D0584K=1 9NSYM 
EG=0.n 

D0533J=1 ,NSYiM 

563 EQ=EO+ROOT( J)*U(K,J) 

564 PRINT1017,K9EQ»REPVEC(K) 

C CALCULATION OF SIGNAL INTENCITIES 
F-RimioOG 
PRINT1026 
D0610 1=1 9 NDIM 
PRIMT1o02,I,H( I ) 

D06lQJ=l,NDIM 

IF( I ,EQ, J) GO TO 610 

SIG=Y{ I,J)*Y{ I »J)*(SIGMA( J9J)-SIGMA( I , I ) +ROOT ( J ) -ROOT ( I ) ) 
SI6 INt=STG/RABAB( I 9 J) 

freq=h( n-H{ J) 
plot=freq*chart 

PRrNTl003, J,FREQ,SIG9SIGINT9PL0T 
610 continue 
9^0 continue 

PRINT1006 
PRlNTlnOb 
GO TO 1 
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ICO:.' F( ?-.'AT(/3H A=,i-10.5»4ri B=,F10«5»4H R=,FlO,5,5H VA=:,F7«3»5H VB = 

1»F7«3//) 

led Format (2H I ,6X,4H£{ n , 8 X» 2 H J,6X,4FSFREQ,9X»3HINT,8X,10HTRANS (I + » 
l»4x» liHTRANS COEFF,* CHART*/) 

iOOC F0RF)AT (I2,F12 o5) 

1003 FORHAT { 20X , 1 2 » 5f 14. 5 3 
lOOA FORMAT ( 2f1&, 8 ) 

1005 FORMAT {* 

1006 F0RNAT(/3 
1007 . FORMAT(8f10,5 ) 

1008 FORMAT( 10F12.5 ) 

1009 FORMAT ( 5 I5^ 

1010 FORMAT ( * MO OF CMS. PER C.P.S. =*,F10.5/ 3 

1011 FORMAT {* DET=*,e12.5/ 3 

1012 FORMATVSH DEH(»I1,4H) = ,El2.5»lOH R00Tr,Il»4H) = ,El2<.5) 

1015 F0R!VAT{4(F18<,8d- = R{*,4ll,^U*) ) 

1016 F0R?lATt29H CALCULATED VALUE OF EQUAT iCN »5X , 15H SUPPLIED VALUE) 

1017 F0RMAT(5X,1H{ ,I2,3H) = » F 1 0 . 6 , 1 5X , F 10 .6 5 

1018 F0Rr'La{40X,*A( I »J)*»50X,*V( I )*/) 
lOiy F0RHAT(8F12.5,6X,F12.53 

1C20 FORMAT! //35X>*SIGMA ZERO MATRIX*/) 

1021 FORMAT {* EQUATION REPLACED=* » I 2 » * DET = * , E 12 . 5 / ) 

1022 FORMAT! 1GX»* MATRIX ELEMENTS OF 1!+)*) 

102:. F0RHAT(2I5,*=**F10.5) 

1024 FORMAT!* C1=*,F10.5,* C2=*,F10.5>* C3=* ,F10 . 5 ) 

1026 FORMAT!* I * ,6X » *E ! I ) * » 8X » * J* , 8X » *FREO* , I IX ,* I NT* » 6X , *INT /R I J I J* , 
18X»*CHART*/ ) 

1030 FORMAT!2i5,F13.8) 

1031 FORMAT!* AS ( * , 1 1 »* »*■* 1 1 » * »* » I 2 » * ) =* ,F 13 . 8 ) 

end 

SIBFTC SUBl NODECK 

subroutine RFDCAL! I»J»K»L»RFD) 

C -HIS SUBROUTINE CALCULATES THE RELAXATION MATRIX ELEMENTS Ro 
C IARTIAL CORsTELATION with either EXTfSSi.iES OF COMPLETE OR NO COREEL AT 1 01 
C IS possible. THE INPUT IS FROM MAIM PROGRAM AI’D THE RESULT IS 

c stored in Red. 

DIMENSION, A! 8 *3 » 15) »F! 15 ) 

COMMON NOSPIN»MQMAX,F 
COMf'ON/ANIL/ A » C1»C2,C3 
N1 IM = 2**NGSPIN 
^iOOP = NQMAX/NOSPIN 
RFD=0.0 

DO 230 NQ=1,NQMAX 

iMQ=NO+NOOP 

IQ=:MQ+.N00P 

LQ=NQ-N00P 

JG=L0-N00P 

RQ=OoO 

RO=RQ+A! I »K>NQ)*A!J»l»NQ)*F!MQ) 

IF ! NQ.LE. NOOP ) RQ=RQ+C1*SQRT ( F ! NQ ) *F ! MQ ) ) *A{ I »fC»NQ )*A ! J»L »M0) + 

lC2*SQRT{F!r'Q)*F! IQ) )*A! I rK»NQ)*A! J»L,IQ> 

IF ( NQ.GT. NOOP .AND.LQoLEoNOOP )RQ=RQ+C3*SQRT ! F ! NQ ) *F !MQ ) ) *A ! I » K.»NQ } * 



Al-10 


1A( J»L,HQ)+Ci-»SURT(F(LO)'“-F(MQ) )-»A( I » K ,;MQ ) «-A ( J » L » LQ ) 

I F ( JQ » GT , NOOP ) RQ = RQ+C2 *A ( I , K , NO ) -“‘A ( J ,L , JG ) ^‘SGR T ( F ( NO ) *F ( JO 5 ) +C3-»- 
loQRT(F(r‘iQ)*F(LO) )-“-A( I »K,NCJ<-A{J,L»LO) 

IF( J-L)220,210»220 
2x0 IF{NQ-NC0P)212»212»211 
211 rF(LO-NOOP) 214,214,216 
2" 2 00213 ■■1=1, NO I'M 

RO = RQ-A ( r ,K .00 ) *A ( M > I » wQ ) «-F ( 00 ) /2 .0 
2 1 3 ivG = OQ-C 1 * A { M , K . ■ iQ ) * A ( M , I , iM-Q ) * SC R T { F ( NQ ) * F ( MQ ) ) / 2 . 0-C 2 *SGR T ( F { NO ) 

IF { IQ ) r1.K,'>iQ )^s-A{M, I , IQ) /2 .0 

GO TO 220 
2x4 D0215 m=1 ,NDIM 

RO = RQ-A(M,K.NQ)-ifA(H. I .WO)*F{riQ) /2.0 

215 Kt =RQ-C3*A (M »K. »NQ ) *4 ( M » I .HQ ) *SGRT ( F ( NO ) *F ( HO ) ) /2 . 0 -C 1*SQRT ( F ( NO) * 
IF(LQ) J*ACM»K.r4Q)-K-A(H>I .LQ)/2,0 

GO TO 220 

216 D02i7H=l ,NDIM 

RQ = RO-A ( M , K . NO ) *A C M » I . NQ )-^F{ NQ ) / 2 . 0 

217 HQ = RG-C2 *A ( H , K . NO ) *A ( H , I , JQ ) --'SOR T ( F ( MO ) *F ( JQ ) ) / 2 . 0-C3*SQR T ( F { NQ ) * 
1F(L0) i*A(M*A»MQ)-''^A(M> I .LQ5/2oO 

220 IF( I-K)230.221 ,230 

221 IF(NQ-NQ0P)223,223,222 

222 IF{LQ-fvOOP)225,225,227 

223 D0224M=1 ,NDIH 

RQ=RO-A(M.t .NQ)*A(M»J»NQ)*F(NQ)/2,0 

224 RQ=RG-C1 * A { 1 . L , NG ) -Jf A ( . i , J , M3 ) -"‘SOR T ( F ( NQ ) *F ( ‘ IQ ) ) /2 « 0-C2*SQR T ( F { NO ) * 
IF ( I C ) ) * A ( , L . NQ ) *A ( M , J , I Q ) / 2 « 0 

GO TO 23o 

225 D0225H=1 ,NDIM 

RQ=nQ-A { M 9 L . NQ ) *A ( M . J . NQ ) *F ( NQ ) /2 „ 0 

226 F,0 = R0-C3-i^A(M,L»N.Q)^A(M, Jj.v.Q) "'^SQRT (F(NG)*F( HQ) ) /2.0-Cl*SQRT { FCNQ)* 
IFIlQ) )*A(M,L»NQ)*A{H9J.LQ)/2oO 

GO TO 230 

227 00228,M=1 ,NDIM 

RQ=RO-A ( H 9 L 9 NO ) ( M 9 J , MG ) *F ( NQ ) / 2 . 0 

228 RQ = RQ-C2* A ( H 9 L 9 NO ) ^>A ( M 9 J 9 JQ ) *SQRT { F ( NQ ) *F ( JQ ) ) / 2 . 0-C3*SQRT ( F ( NQ ) * 
IF ( LO ) 5 ^s-A ( H 9 L 9 NQ ) *A ( M 9 J 9 LQ ) /2 o 0 

■ 23o RFD=RFD+RQ 
RETURN 
END 

‘TIBFTC SUU2 NODECK 

SUBROUTINE DETER(B 9 DET) 

C THIS subroutine CALCULATES THE VALUES OF A DETERMINANT BY PIVOTAL 
C method. no is the order of the DETErF'-UNANT TO be calculated And 
C DET IS its calculated value. 

COM:-ON/KUf*AR/MO 

DlMEMSlOf B( IO 9 IO) 9 A( 10,10) 9PR0W( 10) 

riUM=NO 

DET=l, 

003 1 = 1, NO 
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DC3J=1,mo 
- * ( I » J)=B( I ,J) 

+ J06 = 

^ A ( 1 

CONTinu<^ 

PET=Oo 
GO TO 13 
9 PlVOT=/\{ 

007 I = 

PROW( I ) =/■ ( I ,.Vj j 

/Pivot 

DO 13 I=2,rl(jf4 
DO 13 J=l,NUiv; 

1-'^ I »J5=A( I ,j}-pR0W( I )^Af 3 I 
D020I=2,NUM 
D02uJ=l , nuM 
IF( J-M)13 ,20,10 

Id A( I-l,J-i}=:^( j 
GO TO 20 

A ( I-l 5 J) =A( I , J J 

20 CONTINUE 

num=num-i 

Dp=P I VO^.'fOET^- ( ( -1 . } { M+i , . 

1A 18,16,4 

19 END 


>1 ) }*DET 
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APPEtIDIX II 


H.M.Pl, spectra of 1, 3-Difl^ioro--4 jS-dinitrobenzene 
and SolTent Effects on Coupling Constants. 


Eeprintod froa Eolecular Physics, 1967 
Yol. 12, IIo. 6, 593-596. 
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Reprinted from: 

MOLECULAR PHYSICS, 196?, YOL% 12, Ho. 6, 593-596 

I.M.R. SPECTRA OP 1, 3-DIPLUORO-4,6-DIHITROBEJIZENE AUD 
SOLVENT EPPECTS OH COUPLING CONSTAI'ITS 

by AITIL KUMAR 

Department of Physics, Indian Institute of Technology, 

Kanpur, India 

(Recieved 13 Pebruary 1967) 

1. INTRODUCTION 

The high resolution nuclear magnetic resonance 

spectra and the solvent and concentration dependence of the 

spin-spin coupling constant of 1,3-difluoro 4,6-dinitroben- 

z;ene are reported here. This molecule is a four spin system 

of the type ABXg (iii "tb-e notation of [1]) with nuclei 

labelled as shown in figure 1. 

(B) 

H 

(X)p/2^\P(X) 

(A) 

Figure 1. l,3-difluoro-4-,6~dinitrobenzene. 

Analytical expressions can be derived for all the allowed 
transitions and their relative intensities [2], The 
and ^^P resonance spectra of the molecule, recorded on a 
Varian HR-100 spectrometer operating at 100 B/Ic/s and 94,1 
Mc/s respectively are shown in figure 2. The observed 
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. a 1)10 1. on, A ;,:uid .B, spocbrn of 

l,3-difluoro-4,6-dinitrobc.«..no. Tho 

in. 1 rooonnnoc arc broader oomparod to B 
reounanoo, probably duo to a small coupling 

Of the- A proton with ni,:.rogon, i.,. tho no,-.rl,,Y 
MOg group, whlo.H is rel.TxinH rapidly duo to 
lu.adrupolo intoraotion. x'ho ooparation bet- 
ween the outer linos of tho trlpolots in A 

and B rGscn,-;mcGG is 2j” and Pj“ rnw,,v,f.4-,- ,r 1 

J.3J1 < lospuotivoly, 

In Z resonance the sopar.ation betwoon tbo 

outor linos is rlj® I . | h a 1, 


hetwoon tho inner linos is (|j° 


HP 
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spectra, however, do not show any further structure, 
indicating that the para H-H coupling is less than 0.3 c.p.s. 
('-'half the linewidth) , Thus if is assumed to he zero 

for this system the observed proton and fluorine spectra give 
in a straightforv/ard manner the values of the meta and 
ortho H-F coupling constants, 

2. SOLVEIf EFIBCIS 

1 , 3-difluoro-4- , 6~dinitrobenzene is a solid at 
room temperature and can he studied in high resolution ■ 
through the use of a solvent, Ihe spectra were studied in 
various solvents to see the effect, if any, on the spin 
coupling constants due to the solvent. The results are 
given in the table. In two of the solvents deutero- 
chloroform and dioxane the spectra were studied at several 
different concentrations. The variation of the coupling 
constants vrith concentration, in the range of concentrations 
studied, is smaller than the variation in different solvents. 
The coupling constants tend to increase Yri.th percentage 
concentration in d-chloroform. This trend is reversed in 
dioxane and further the values are smaller in d-chloroform 
than in dioxane. The coupling constants in the pure 
compound therefore are expected to lie between the values 
in these two solvents. 

The solvent and concentration dependence of 
coupling constants, except where they are due to configura- 
tional modifications, indicate changes in intramolecular 




standard deviation for. about ten or more readings. 


Dieleotric constants from; 1961, Handbook of Ohemlstry and Physics (Ohio; 
Chemical Rubber Co.), fkat of THP from: Meloan, C, A., and Kiser, R.W. , 
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interactions. The Tarious factors considered responsible 
for these changes in and coupling 

constants, reported by several workers [3-8] » include, the 
reaction field effects of the medium, molecular association, 
complex formation and the shape of the solvent and solute 
molecules. The reaction field [9] is an electric field 
produced at the solute molecule due to a polarization of . 
the medium by the polar solute. The effect of this is 
indicated by similar changes in chemical shifts and coupling 
constants with dielectric constant of the solvent [5-8], 


The present observed variation (table) of the sum 
of the two coupling constants at 10 per 

cent concentration in the first four solvents is approxi-^ 
mately linear with dielectric constant. This variation is 
therefore likely to be due to the reaction field effects of 
the medium [3-8], In the other three solvents namely 
acetonitrile, tetraliydrofuran and dioxane, the coupling 
constants bear no such relation to the dielectric constants. 
These three later solvents are known to be high donor centres 
[10] and will probably form donor-acceptor complexes by 


-j- This s-um, which can also be obtained directly from the 

fluorine spectra, shows a larger variation between solvents 
than either or separately. The similar 


HF 

variations of J 


HF 

° and IxT-m™ with solvents indicate that 


these two coupling constants are of the same sign [6]. 
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hydrogen bonding between the ring protons of the solute 
molecule and the donor centres of the solrent molecules 
[4,7]. The complex formation may affect the coupling 
constants in a way which is different from the reaction 
field. 

Both the reaction field and the complex formation 
are also knom to affectthe chemical shifts [8,11] and 
therefore it is reasonable to expect some correlation 
between the variation of the coupling constants and that 
of the chemical shifts.. An approximate correlation is 
obtained in the present case. The two coupling constants 
tend to increase with increasing chemical shifts, and 
~\/-Q and decreasing values of ( 1i^) in agreement with 

previous observations [3-6]. The chemical shifts in 
benzene are, however, very different from those in other 
solvents, - probably because benzene is an aromatic compound. 

These observations suggest that the effect of the 
solvent on the coupling constants should be taken into 
account, particularly ^ in solid samples where solvents are 
unavoidable.. The apparent concentration independence of 
coupling constants within experimental errors could, in . 
several cases, be due to the limited range of concentrations 
which could be studied (limited by solubility, signal-to- 
noise ratio, etc.) rather than the absence of any solvent 
effects. Therefore, the solvent dependence of coupling 
constants should be determined by studying the concentration 
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dependence in at least two solvents of different natures, 
and the pure compound values obtained by extrapolation, 
just as is the practice for chemical shifts. Conclusions 
arrived at otherwise could be erroneous [12], The study 
of solvent dependence of the coupling constants should also 
throw some light on the Various intramolecular interactions 
responsible for the coupling between two spins in a ' 
molecule, 

She author wishes to acknowledge many invaluable 
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interest taken by Professor P. ¥enkateswarlu in this work. 
Thanks are also dxie to hr, Eamesh Ohand for sending the 
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REPEEEHCES 

[1] Pople, J.A., Schneider, W.G. , and Bernstein, H,J,, 
1959,'*High Resolution ETuclear Magnetic Resonance]^ 
Chap. 6 (Hew York; McGraw-Hill), 

[2] Mortimer, E,S,, 1959» J, Mol. Spectrosc., _3,. 335. 
Cohen, A.h,, and Sheppard, H. , 1959, Proc. R. Soc, 
k,, 252 , 488, de Kowalewskl, h.G, ,. and Kowalewski, 
V.J., 1961, J. Phys. Radium, Paris, 129, 

[3] Watts, V.S,, and Goldstein, J.H., 1963, J, Mol, 
Spectrosc., U, 325; 1965, J, Chem, Phys,, 4_2, 228. 
Shapiro, B.L,, Ebersole, S,J., and Kopshik, R.M. , 
1963, J. Mol. Spectrosc., U,, 326 ; 1963, J. Chem. 
PHys., 29, 3154. 



AII-ll 


[4] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 


[ 11 ] 

[ 12 ] 


oLT’p!" S.S., 1964, J. . 

OJiem. Phys., 2415 . » d. 

Hutton, H.M., and Schaefer, f, 1965 Cari t n-h 
12 . 3116. Fg, s T ^ 

1961 r nJ ®’ ^•’ Sederholm, O.H., 

1964, 1 . cheM. ays., 2090; 1965, iba.. 42, 79 . 

Smith, S.I., and Cox, E.H., 1965 j MoT + 

16, 216- T n-u ^7 Spectrosc., 

215, 1966, J. chem. Phys., 45, 2848. 

■ TrEmkiss, S.S., 1963, j. ^ 

Bvatis, D.P., 1963 , J. ohem. Soc., p. 5575 . 

Schneider, W.G., 

W60,^J. Chen. Phys., 1227 , i 960 . Can. l.^ben., 

honplene (Berlin, Gottingen, Heidelberg: Sprier- 

Reference [l], chap. I 5 . 

Hageswara Rao, B.P 

„ , ^ j-yb 5 , Molec, Pnys,- 7 . 507 

Schaefer, 1 . , 1962 Cm t nv ”* 

t , xi?o,d, oan. J, Chem,, 431. 



